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Abstract 



We review spacetime dynamics in the presence of large-scale electromagnetic fields 
and then consider the effects of the magnetic component on perturbations to a spa- 
tially homogeneous and isotropic universe. Using covariant techniques, we refine and 
extend earlier work and provide the magnetohydro dynamic equations that describe 
inhomogeneous magnetic cosmologies in full general relativity. Specialising this sys- 
tem to perturbed Friedmann-Robertson- Walker models, we examine the effects of 
the field on the expansion dynamics and on the growth of density inhomogeneities, 
including non-adiabatic modes. We look at scalar perturbations and obtain analytic 
solutions for their linear evolution in the radiation, dust and inflationary eras. In 
the dust case we also calculate the magnetic analogue of the Jeans length. We then 
consider the evolution of vector perturbations and find that the magnetic presence 
generally reduces the decay rate of these distortions. Finally, we examine the im- 
plications of magnetic fields for the evolution of cosmological gravitational waves. 

PACS: 98.80.-k; 98.62.En; 98.65.Dx 
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1 Introduction 



Magnetic fields are a widespread and significant component of the Universe [I]-[7j. The Milky 
Way and many other spiral galaxies possess magnetic fields with strengths of order a few 
micro Gauss, that are coherent over the plane of the galactic disc. Magnetic fields are also a 
common property of galaxy clusters, extending well beyond their core regions. The strengths 
of the ordered magnetic fields in the intracluster medium exceed those typically associated 
with the interstellar medium of the Milky Way, suggesting that galaxy formation, and even 
cluster dynamics, could be influenced by magnetic forces. Furthermore, reports of Faraday 
rotation in high-redshift, Lyman-a absorption systems suggest that dynamically significant 
magnetic fields could be present in protogalactic condensations. 

Despite their widespread presence, however, the precise origin of any large-scale cosmological 
magnetic field is still a mystery, and its likelihood or necessity remains a subject of debate [5]- 
[12]. The alignment of the galactic fields, especially those seen in spiral galaxies, seems to 
support the dynamo idea [13]- [T7]. However, the galactic dynamo cannot work without an 
initial seed field [18j-|20j and these seeds may require I?-fields of primeval origin [21]-[H]. De- 
ciding whether galactic and cluster magnetic fields are primordial relics or post-recombination 
artefacts is difficult because the strong amplification of the fields in these virialized systems 
has overwhelmed all traces of their earlier history. In contrast, possible magnetic imprints 
in the Cosmic Microwave Background (CMB), or a magnetic field in the intercluster space, 
can provide better insight into the early history because the presence of small-scale magnetic 
fields leaves small angular scale features in the CMB undamped [HUB], and leads to distinc- 
tive polarisation signatures |47| . The idea of primordial magnetism is attractive because it can 
potentially explain all the large-scale 5-fields seen in the universe today, especially those found 
in high-redshift protogalaxies. Early magnetogenesis is not problem-free however. Seed fields 
generated during, say, phase-transitions in the radiation era have typical coherence lengths 
that are very small and will destabilize the dynamo. On the other hand, fields which survived 
an epoch of inflation are typically too weak to sustain the dynamo. Although the literature 
contains several mechanisms of primordial magnetic amplification [IB]- [21], with early reheat- 
ing looking like a promising stage [58j, the issue remains open. Any cosmological magnetic 
field must also be consistent with a number of astrophysical constraints. These come from 
primordial nucleosynthesis [59 60J, and the high isotropy of the CMB radiation [61]-[80]. This 
constrains the current strength of a large-scale homogeneous 5-field below ~ 10~ 9 G [BB] (see 
also § 111.41 here), although tangled random fields on small scales could reach up to ~ 1CT 6 G. 

Large-scale magnetic fields of fiG strength can affect the evolution of structure in the Universe, 
and studies of their effects have a long history. Anisotropic, spatially homogeneous, relativis- 
ts cosmologies permeated by large-scale magnetic fields have been analysed in [EE]- [S3] and 
more recently in [B2] and [HI]- [EH]- The bulk of the available inhomogeneous treatments are 
Newtonian, with general relativistic treatments only a recent addition to the literature. The 
evolution of linear density and vorticity perturbations in magnetic Newtonian cosmologies 
were addressed by [90]- [93], and again more recently by [94195] ; general relativistic treatments 
have been given in [96j-|98j. These investigations have established that magnetic fields are 
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sources of density and vorticity distortions, but the rather complicated action of the 5-field 
did not allow for analytic solutions to the relativistic equations. Solutions were provided with 
the help of covariant and gauge-invariant techniques, which considerably simplify the math- 
ematics of cosmological magnetohydrodynamics [5§]- [lU2j . The covariant equations were also 
used as the basis for qualitative studies by means of phase-plane methods |103j . while an 
analogous approach has been used to analyse two-component charged plasmas |104|ll05j . In 
addition to the linear regime, a number of authors have looked into the non-linear interac- 
tion between cosmic magnetic fields and electrically-charged fluids, primarily on sub-horizon 
scales. This includes work on the amplification of the 5-field by shearing effects during the 
anisotropic collapse of a proto-galactic cloud [1U6] - |116] . The aforementioned studies combine 
analytical and numerical approaches, and apply to the pre- or to the post-recombination era of 
the universe. All this has motivated and facilitated work on the potential effects of magnetism 
on the CMB spectrum, including its polarisation [THl63] -[80]. and the low-quadrupole moment 
problem |117j . 

In the present article, we use covariant and gauge-invariant techniques to analyse the effects of 
cosmological magnetic fields on the large-scale structure of the universe. We do so by refining 
and extending the work of [99]-p.01J. After a brief introduction to the 1+3 covariant formalism 
in §[2J we use §[3], [Hand [5] to provide a detailed discussion of cosmological electrodynamics and 
magnetohydrodynamics in full general relativity. Section [H] utilises this information to study 
the effects of the 5-field on weakly magnetic almost-Friedmann-Robertson- Walker (FRW) 
cosmologies in a gauge-invariant manner. One of the key results is that linear perturbations in 
the magnetic energy density evolve in step with those in the density of the matter. Focusing 
on scalar perturbations, we look at the evolution of linearized matter perturbations during 
the radiation and dust eras in § During the radiation era, and on super-horizon scales, 
our solutions show the magnetic pressure slows down the standard 'Jeans' growth rate of the 
density contrast. At the same time, the 5-field increases the oscillation frequency of small- 
scale density perturbations in proportion to its strength relative to the radiation background 
density. The same amount of increase is also observed in the associated Jeans length. After 
the radiation era ends, the fluid pressure is effectively zero, and the magnetic field emerges as 
the sole source of pressure support. We calculate the associated 'magnetic' Jeans length and 
find that for fiG fields its scale is comparable to the size of a galaxy cluster. Given that during 
the dust era the magnetic relative strength decays with time, we also calculate the late-time 
magnetic effects on density perturbations. Here too, we find that the magnetic pressure in- 
hibits the gravitational clumping of matter. This section closes by looking at the evolution 
of magnetic density perturbations within an inflationary false- vacuum environment. We find 
that the outcome is very sensitive to the effective equation of state of the cosmic medium. 
Section [8] looks at the magnetic effects on isocurvature density fluctuations, namely on pertur- 
bations evolving on irrotational spatial hypersurfaces of constant curvature. After establishing 
the necessary conditions for the existence of such disturbances in the presence of the .B-field, 
we find that they always decay in time. In § [HI we treat the magnetic medium as a two-fluid 
system and consider the effective entropy perturbations coming from differences in the dy- 
namical behaviour of the two components. Our study shows that such disturbances vanish 
during the radiation era, while after recombination they either vanish, or remain constant, 
depending on the initial conditions. The magnetic effects on rotational (vector) density per- 
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turbations are considered in § [TUJ At the linear perturbative level these distortions are directly 
related to the amount of vorticity in the universe. Focusing on the dust era, we show that 
the magnetic presence reduces the standard decay rate of rotational perturbations. This leads 
to more residual vorticity relative to magnetic-free universes. The implications of the 5-field 
for the evolution of gravitational waves are examined in § [TTJ As a first step, we introduce 
the additional constraints needed to isolate the linear pure-tensor modes in the presence of a 
magnetic field. Once these have been imposed, we consider the magnetic effects on the propa- 
gation of gravitational-wave distortions. We also look at the role of the magnetic anisotropic 
pressure during the radiation era and discuss the critical, zero eigenvalue, case that emerges 
there. 



2 1+3 covariant description 



The covariant approach to general relativity and cosmology dates back to the work of Heck- 
mann, Schucking, and Raychaudhuri in the 1950s [118|119] and it has since been employed 
in numerous applications by many authors. The formalism uses the kinematic quantities of 
the fluid, its energy density and pressure and the gravito-electromagnetic tensors instead of 
the metric, which in itself does not provide a covariant description. The key equations are the 
Ricci and Bianchi identities, applied to the fluid 4-velocity vector, while Einstein's equations 
are incorporated via algebraic relations between the Ricci and the energy-momentum tensors. 
Here, we will only give a brief description of the approach and direct the reader to a number 
of review articles for further details and references [120] - [T23] . 



2. 1 Local spacetime splitting 



Consider a general spacetime with a Lorentzian metric g a b of signature (— , +, +, +). Introduce 
a family of fundamental observers along a timelike congruence of worldlines tangent to the 
4-velocity vectoiQ 

u a = ^, (2.1.1) 

where r is the associated proper time and u a u a = —1 [122] . This fundamental velocity field 
introduces a local, 1+3 'threading' of the spacetime into time and space. The vector u a de- 
termines the time direction, and the tensor h ab = g ab + u a Ub projects orthogonal to 4-velocity 
field into what is known as the observer's instantaneous rest space. In the absence of rotation, 
u a is hypersurface orthogonal and h ab acts as the metric of the 3-D spatial sections. 

Employing u a and h abl we define the covariant time derivative and the orthogonally projected 



1 Latin indices run between and 3 and Greek vary from 1 to 3. We use geometrized units with 
c = 1 = 8ttG, which means that all geometrical variables have physical dimensions that are integer 
powers of length. 
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gradient of any given tensor field S ab .. according to 

S ab ... cd - = a'Y.S,,,,..:-' and D e S ab .. cd - = h e s hj h b p h q c h r d ■ ■ ■ V.,.^,...'"'" , (2.1.2) 

respectively. The former indicates differentiation along the timelike direction and the latter 
operates on the observer's rest space. 

2.2 Matter fields 

Relative to the fundamental observers, the energy-momentum tensor of a general imperfect 
fluid decomposes into its irreducible parts as [122J 

Tab = pUaUb + PKb + 2q(aUb) + tab • (2-2.1) 

Here, p = T ab u a u b and p = T ab h ab /3 are, respectively, the energy density and the isotropic 
pressure of the medium, q a = —h b T bc u c is the energy-flux vector relative to u a , and 7i a b = 
h( a c hb) d Tpiis the symmetric and trace-free tensor that describes the anisotropic pressure of 
the ffuidlfj It follows that q a u a = = ir ab u a . When the fluid is perfect, both q a and ir ab are 
identically zero, and the remaining degrees of freedom are determined by the equation of state. 
For a barotropic medium the latter reduces to p = p (p), with c 2 s = dp/ dp giving the the square 
of the associated adiabatic sound speed. 

When dealing with a mult i- component medium, or when allowing for peculiar velocities, one 
needs to account for the velocity 'tilt' between the matter components and the fundamental 
observers (e.g. see [124] - [T28] ). Here we will consider a single-component fluid and we will 
assume that the fundamental observers move with it. 

2.3 Covariant kinematics 

The observers' motion is characterized by the irreducible kinematical quantities of the Ma- 
congruence, which emerge from the covariant decomposition of the 4-velocity gradient 

V b u a = a ab + u ab + ^ Qh ab - A a u b , (2.3.1) 

where a ab = D/ b u a \, u) ab = D[ 6 w a i, 6 = V a u a = D a u a and A a = u b V b u a are respectively 
the shear and the vorticity tensors, the volume expansion (or contraction) scalar, and the 4- 
acceleration vector p_22j . Then, a ab u a = = u) ab u a = A a u a by construction. The volume scalar 
is used to introduce a representative length scale (the cosmological scale factor a) by means of 
definition a/a = 6/3. Also, on using the orthogonally projected alternating tensor e abc (with 

2 Angled brackets denote the symmetric and trace-free part of second-rank tensors projected or- 
thogonally to u a and the projected component of vectors (i.e. = h a b v b ). 
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£ a bc = 3u[ a ebc]dA d - see [135J ) , one defines the vorticity vector u a = e a b c uj bc /2. Note that e a b c = 
r} a bcd,u d , where r\ abcd is the totally antisymmetric permutation tensor of the spacetime. This is a 
covariantly constant quantity, with r s] a b c sf^' p ' i = — A\5y a e 8j Sj'djf and 77 0123 = [— &et(g ab )Y l l 2 . 
The tensor v a b = D fe w a = a a b + u> ab + (Q/S)h ab describes the relative motion of neighbouring 
observers (with the same 4- velocity). In particular, v a = v a bX b is the relative velocity of the 
associated worldlines and \a is their relative position vector (see [121|I129] for details). 

The non-linear covariant kinematics is determined by a set of three propagation equations 
complemented by an equal number of constraints [122] . The former contains Raychaudhuri's 
formula 

@ = ~l ^ " \ (P + 3P) ~ 2((j2 " ^ + DMfl + ^ ' (2 - 3 - 2) 
for the time evolution of G; the shear propagation equation 

2 c 1 

<j{ab) = — g ©^afe - 0"c{a0" 6 b) ~ ^(a^b) + D{aA> + A( a A b ) - E ab + - 7l ab , (2.3.3) 

which describes kinematical anisotropics; and the evolution equation of the vorticity 

2 1 

U( a ) = -- - - curl A a + a ab u b . (2.3.4) 

Here a 2 = cr a bcr ab /2 and u 2 = u a bU} ab /2 = u a u a are respectively the scalar magnitudes of 
the shear and the vorticity, while E ab is the electric component of the Weyl tensor (see § 14.21 
below). Also, cxulva = e a bcD b v c for any orthogonally projected vector v a . 

Equations (12.3.21) . (12.3.31) and (I2.3.4[) are complemented by a set of three non-linear constraints. 
These are the shear 

D b a ab = - D a + curl u a + 2e abc A b u c - q a , (2.3.5) 

the vorticity 

D a uj a = A a u a , (2.3.6) 

and the magnetic Weyl constraint 

H ab = cur\a ab + B {a uj b) + 2A {a u b) , (2.3.7) 

where curlT a fc = e c d(c£) c T b ) d for any symmetric, orthogonally projected tensor T ab . We finally 
note that the cosmological constant has been set to zero (i.e. A = 0) throughout this review. 



3 Electromagnetic fields 

Covariant studies of electromagnetic fields date back to the work of Ehlers [12U] and Ellis [129J. 
In addition to its inherent mathematical compactness and clarity, the formalism facilitates a 
physically intuitive fluid description of the Maxwell field. This is represented as an imperfect 
fluid with properties specified by its electric and magnetic components. For a fully covariant 
study of electromagnetic fields in curved spacetimes the reader is referred to [130] . 



S 



3.1 Electric and magnetic components 



The Maxwell field is covariantly characterized by the antisymmetric electromagnetic (Faraday) 
tensor F ab . Relative to a fundamental observer, the latter decomposes as [129] 

E ab = 2u [a E b] + e abc B c , (3.1.1) 

where E a = F ab u b and B a = e abc E bc /2 are respectively the electric and magnetic fields mea- 
sured by the observer. Note that E a u a = = B a u a , ensuring that both E a and B a are spacelike 
vectors. Also, B a = e abc E bc /2 guarantees that B a is the dual of F ab . 

The Faraday tensor also determines the energy-momentum tensor of the Maxwell field accord- 
ing to 

T^ m) = ~F ac F c b - 1 F cd F cd g ab . (3.1.2) 

The above combines with (I3.1.ip to facilitate the irreducible decomposition of T a ^ m ^ relative 
to the u a - frame |129j . 

T^ m) = \{E 2 + B 2 ) u a u b + \(E 2 + B 2 ) h ab + 2Q (a u b) + V ab . (3.1.3) 

Here E 2 = E a E a and B 2 = B a B a are the magnitudes of the two fields, Q a = e abc E b B c is the 
electromagnetic Poynting vector and V ab is a symmetric, trace-free tensor given by 

V ab = V {ab) = \{E 2 + B 2 ) h ab - E a E b - B a B b . (3.1.4) 

Expression (13.1.31) provides a fluid description of the Maxwell field and manifests its generically 
anisotropic nature. In particular, the electromagnetic field corresponds to an imperfect fluid 
with energy density (E 2 + B 2 )/2, isotropic pressure (E 2 + B 2 )/6, anisotropic stresses given by 
Vab and an energy-flux vector represented by Q a - Equation (13.1.31) also ensures that Tjf m ^ a = 0, 
in agreement with the trace-free nature of the radiation stress-energy tensor. Finally, we note 
that by putting the isotropic and anisotropic pressure together one arrives at the familiar 
Maxwell tensor, which assumes the covariant form 

M ab = \{E 2 + B 2 ) h ab - E a E b - B a B b . (3.1.5) 

3.2 Maxwell's equations 

We follow the evolution of the electromagnetic field by means of Maxwell's equations. In their 
standard tensor form these read 

V [c F ab] =0 and V b F ab = J a , (3.2.1) 
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where (13.2. lb ) reflects the existence of a 4-potential and J a is the 4-current that sources 
the electromagnetic field. With respect to the w a -congruence, the 4-current splits into its 
irreducible parts according to 

Ja = PeU a + Ja, (3.2.2) 

with p e = —J a u a representing the measurable charge density and J a = h a h Jb the orthogonally 
projected current (i.e. J a u a = 0). 

Relative to a fundamental observer, each one of Maxwell's equations decomposes into a time- 
like and a spacelike component. Projecting (13.2. la ) and (13.2.1b ) along and orthogonal to the 
4-velocity vector u a , one obtains a set of two propagation equations [129] 



E {a) = (a ab + e ahc u c - - Qh ab ) E b + e abc A b B c + curl B a - J a , (3.2.3) 

B {a) = (a ab + e abc u c - - Qh ab ^j B b - e abc A b E c - curl E a , (3.2.4) 

and the following pair of constraints 

D a E a = p e -2u a B a , (3.2.5) 

D a B a = 2u a E a . (3.2.6) 



Note that, in addition to the usual 'curl' and 'divergence' terms, Eqs. (I3.2.3l) - ()3.2.6p also 
contain effects triggered by the relative motion of neighbouring observers (with the same 4- 
velocity - see § 12.31) . These are carried by the kinematic terms in the right-hand side of the 
above. Thus, p e = —2u a B a is an effective electric charge caused by the relative motion of the 
magnetic field, while 2uj a E a acts as an effective magnetic charge triggered by the relatively 
moving E-G.e\d. The acceleration terms in (I3.2.3P and (13.2.41) . on the other hand, also reflect 
the fact that spacetime is treated as a single entity. 

3.3 Conservation laws 

The twice contracted Bianchi identities guarantee the conservation of the total energy mo- 
mentum tensor, namely that V b T ab = 0. This constraint splits into a timelike and a spacelike 
part, which respectively lead to the energy density and the momentum-density conserva- 
tion laws. The energy momentum tensor of the electromagnetic field satisfies the constraint 
ybjjem) _ _p^jb^ w ith the Faraday tensor given by (13.1.11) and the quantity in the right- 
hand side representing the Lorentz 4-force. Thus, for charged matter the conservation of the 
total energy-momentum tensor T ab = + leads to the formulae 

p = -e(p +p)- D a q a - 2A a q a - a ab 7r ab + E a J a , (3.3.1) 

for the energy density, and 
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(p + p)A a = -D a p - q {a) --9q a - (cr a b + ^ ab )q b - B b n ab - n ab A b 
+p e E a + e abc J b B c , 



(3.3.2) 



for the momentum density. The last two terms in the right-hand side of (I3.3.2P represent the 
familiar form of the Lorentz force. We also note that the electromagnetic effects depend on 
the electrical properties of the medium (see § 13.41 below). 

The antisymmetry of the Faraday tensor (see Eq. (13.1.11) ) and the second of Maxwell's formulae 
(see Eq. (13. 2. lb )) imply V a J a = and therefore ensure the conservation of the 4-current 
density. Using decomposition (I3.2.2p . we arrive at the covariant form of the charge- density 
conservation law [129] 

p e = -Q Pe -D a J a -A a J a . (3.3.3) 
Hence, in the absence of spatial currents, the evolution of the charge density depends entirely 
on the volume expansion (or contraction) of the fluid element. 

3.4 Ohm's law 

The relation between the 4-current and the electric field, as measured by the fundamental 
observers, is determined by Ohm's law. Following [1311)132] . the latter has the covariant form 

Ja = PeUa + SE a , (3.4.1) 

where <; is the scalar conductivity of the medium. Thus, Ohm's law splits the 4-current into 
a timelike convective component and a conducting spacelike counterpart. Projecting (13.4.11) 
into the observer's rest space gives 

J a = ^E a . (3.4.2) 
This form of Ohm's law covariantly describes the resistive magnetohydrodynamic (MHD) 
approximation in the single-fluid approach. Note the absence of the induced electric field from 
the above, reflecting the fact that the covariant form of Maxwell's formulae (see expressions 
f l3.2.3l) - fl3.2.6p ) already incorporates the effects of relative motion. According to (I3.4.2p . non- 
zero spatial currents are compatible with a vanishing electric field as long as the conductivity of 
the medium is infinite (i.e. for q — > oo). Thus, at the limit of ideal magnetohydrodynamics, the 
electric field vanishes in the frame of the fluid. On the other hand, zero electrical conductivity 
implies that the spatial currents vanish, even when the electric field is non-zero. The electrical 
conductivity is typically treated in a phenomenological manner and here we will also assume 
that it remains constant throughout the medium. 



4 Gravitational field 

Covariantly, the local gravitational field is described by a set of algebraic relations between the 
Ricci curvature tensor and the energy-momentum tensor of the matter. The free gravitational 
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field, on the other hand, is described by the electric and magnetic components of the conformal 
curvature (Weyl) tensor. 



4-1 Local Ricci curvature 



In the general-relativistic geometrical interpretation of gravity, matter determines the space- 
time curvature, which in turn dictates the motion of the matter. This interaction is evident 
in the Einstein field equations, which in the absence of a cosmological constant take the form 

Rab = T ab -^ Tg ab , (4.1.1) 

where R ab = R c ac b is the spacetime Ricci tensor, T ab is the energy-momentum tensor of the 
matter fields, with T = T a a being the trace. For our purposes the total energy-momentum 
tensor has the form T ab = + T^ m) , where is given by Eq. (12X11 and T a ^ m) by fl3X3|) . 
Thus, 



T, 



ab ' 



p+- (B' + E< 

-TTab + Vab , 



U a U b + 



P+\(B 2 + E< 



h ab + 2(g (a + Q(a)ub) 



(4.1.2) 



ensuring that p + (B 2 + E 2 )/2 is the total energy density of the system, p + [B 2 + E 2 )/6 is 
the total isotropic pressure, q a + Q a is the total heat-flux vector and 7i a b + V a b is the total 
anisotropic pressure. The inclusion of electromagnetic terms in the energy-momentum tensor 
of the matter guarantees that the contribution of the Maxwell field to the spacetime geometry 
is fully accounted for. 

The successive contraction of the Einstein field equations, assuming that T ab is given by 
Eq. (14.1.21) . leads to the following algebraic relations: 



RabU a U b . 

haR bc u c - 



p 



3p + E 2 + B 2 ) 



(q a + Qa 



h n c h b R, 



cd • 



1 

3 



i>-n+ - ( /: 2 • r>- 



i Lab 



ab 



(4.1.3) 
(4.1.4) 
(4.1.5) 



In addition, the trace of f)4.1.ip gives R = — T, with R = R a a and T = T a a = 3p — p, 
where the latter result is guaranteed by the trace- free nature of T^^. Note that the above 
expressions are valid irrespective of the strength of the electromagnetic components and recall 
that q a = = n ab when dealing with a perfect fluid. 
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4-2 Long-range Weyl curvature 



The Ricci tensor describes the local gravitational field of the nearby matter. The long-range 
gravitational field, namely gravitational waves and tidal forces, propagates through the Weyl 
conformal curvature tensor. The splitting of the gravitational field into its local and non-local 
components is demonstrated in the following decomposition of the Riemann tensor, 

1 1 

Rabcd = C a bcd + ~ {9acRbd + 9bdRac ~ 9bcRad ~ 9adRbc) ~ ~ R \9ac9bd ~ 9ad9bc) , (4-2.1) 

2 o 

where C abcd is the Weyl tensor. This shares all the symmetries of the Riemann tensor and 
is also trace-free (i.e. C c acb = 0). Relative to the fundamental observers, the Weyl tensor 
decomposes into its irreducible parts according to 

C a bcd i^9abqp9cdsr VabqpVcdsr*) U^U E^ ( K T]abqp9cdsr 9abqpVcdsr S ) U^U , (4.2.2) 

where r\ abcd = V[abcd] is the spacetime permutation tensor defined in § 12.31 and g a bcd = 9ac9bd — 
9ad9bc (e.g. see |133|ll34j ). The symmetric and trace- free tensors E ab and H ab are known as the 
electric and magnetic Weyl components and they are given by 

E ab = C acbd u c u d and H ab = ^ e a cd C cdbe u e , (4.2.3) 

with E ab u b = = H ab u h . Given that E ab has a Newtonian counterpart, the electric part of 
the Weyl tensor is associated with the tidal gravitational field. The magnetic component, on 
the other hand, has no Newtonian analogue and is therefore primarily associated with gravi- 
tational waves and spatial 3-curvature anisotropy |122j . Of course, both tensors are required 
if gravitational waves are to exist. 

The Weyl tensor represents the part of the curvature that is not determined locally by matter. 
However, the dynamics of the Weyl field are not entirely arbitrary because the Riemann tensor 
satisfies the Bianchi identities. When contracted, the latter take the form [133J 

V d C abcd = V [b R a]c + - g c[b V a] R , (4.2.4) 

by means of decomposition (I4.2.ip . In one sense the contracted Bianchi identities act as the 
field equations for the Weyl tensor, determining the part of the spacetime curvature that 
depends on the matter distribution at other points [133] . The form of the contracted Bianchi 
identities guarantees that once the electromagnetic contribution to the Ricci curvature has 
been incorporated, through the Einstein field equations, the effect of the Maxwell field on the 
Weyl curvature has also been fully accounted for. 

The 1+3 splitting of (14.2.41) provides a set of two propagation and two constraint equations for 
the evolution of the long-range gravitational field, namely of tidal forces and gravity waves. 
In particular, on using the decomposition (I4.2.2p . the timelike component of (I4.2.4P leads 
to IT221T351 



13 



E 



(ab) 



and 



H 



(ab) 



-QE, 



ab 



p + p + -(B 2 + E 2 



a a b + curl H ab - - (jr ab + V, 



ab 



+3a <a c 



+ 2E cd ( a A c H b ) c 



~'cd{a' 



U) ' 



1 3 
-QH ab - cui\E ab + -(curl7r ab + curlP ab ) + 3a {a c H b)c - - uj( a (q b) 



b), 



~<7 c 6> (g d +g d )-^ c #6 



(4.2.5) 



(4.2.6) 



Taking the time derivatives of the above one arrives to a pair of wavelike equations for the 
electric and the magnetic parts of the Weyl tensor, showing how curvature distortions prop- 
agate in the form of gravitational waves like ripples in the spacetime fabric. These waves 
are also subjected to a set of constraints, which emerge from the spacelike component of the 
decomposed Eq. fl4.2.4j> and are given by [122|I135] 



D b E ab = -B 



p+ l -(b 2 + e 2 



—2>H ab uj b + e abc 



a 



b JJ°d 



D b (Tlab+Vab) 



UJ b (q c +Q c ) 



\®{q a +Qa) + \vab{ 



f + Q b ) 

(4.2.7) 



and 



D b H, 



ab ' 



p + p+-(b 2 + e< 



UJn 



curl (q a + Q a ) + 3E ab u b - - (n ab + V ab ) u b 



—£abc<? d 



E cd + _ (^cd + pcd 



(4.2.8) 



respectively. The above expressions are similar to Maxwell's formulae, which explains the 
names of E ab and H ab . In fact, this Maxwell-like form of the free gravitational field underlines 
the rich correspondence between electromagnetism and general relativity, which has been the 
subject of theoretical debate for many decades (see [T36]- [T39] for a representative list). 



4-3 Spatial curvature 



When the fluid is irrotational, the rest spaces of the fundamental observers mesh together to 
form spacelike surfaces orthogonal to their worldlines. These are normal to the w a -congruence 
and define the hypersurfaces of simultaneity of all the comoving observers. In the presence 
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of vorticity, however, Frobenius' theorem forbids the existence of such integrable surfaces 
(e.g. see |140lll41j ). The w a -congruence is no longer hypersurface orthogonal. Then, one can 
still talk about the observers' spatial rest-space but only locally. The local 3-Riemann tensor, 
defined by 

Tlabcd = h a q h b s hJh d p R qsfp - v ac v bd + v ad v bc , (4.3.1) 
where v ab = ^> b u a is the relative position tensor. On using Eqs. (j4.1.3p - (j4.1.5p and expressions 
fl4XTD . (142^1) . definition fliXTD gives 



1 



P+ l -(E 2 + B 2 )-\& 



(h ac h bd — h ad h bc ) 



+ ~ [haci^bd + V bd ) + {ir ac + V ac )h bd - h ad (7l bc + V bc ) - (tt ad + Vad)h bc } 

~ © [Kc{.cf bd + LO bd ) + (a ac + uj ac )h bd - h ad (a bc + u bc ) - (a ad + LO ad )h bc ] 

-(a ac + uj ac )(a bd + u bd ) + (a ad + w ad )(a 6c + u bc ) , (4.3.2) 

which offers an irreducible decomposition of the local 3-Riemann tensor orthogonal to u a . 
Then, one can easily show that 

Tlabcd = T^[ab][cd] (4.3.3) 

and that 



Habcd — Tlcdab — ~ ~ @ (h ac U bd + UJ a ch b d ~ h ad UJ bc — U ad h bc ) 

—2 {<J ac uj bd + u ac a bd — <J ad u) bc — u) ad a bc ) ■ 



(4.3.4) 



Therefore, in the absence of vorticity, 7Z abcd = lZ cdab and the spatial Riemann tensor possesses 
all the symmetries of its 4-D counterpart. 

In analogy with R ab and R, the local Ricci tensor and Ricci scalar of the 3-D space orthogonal 
to u a are defined by 



/v afe ,L '^cadb /v acb 



and 



n = h ab n ab , 



(4.3.5) 



respectively. Thus, contracting f)4.3.2p along the first and third indices we arrive at what is 
usually referred to as the Gauss-Codacci equation: 



K 



ab 



E a b + ~ 



p + U E 2 + B 2 )-\e 2 + a 2 -uJ 2 



hab + ^ ( n cLb + V a b) 



• - @(o"ab + u ab ) + a c { a a c b ) — uj c ( a u c b ) + 2a c [ a uj c b ] . 



(4.3.6) 



A further contraction leads to the trace of the above and provides the local Ricci scalar of the 
3-dimension space orthogonal to u a 



K = h ab TZ 



ab 



P + \(E 2 + B 2 )-\e 2 + o 



2 u 2 



(4.3.7) 
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This expression is of major importance, since it is nothing else but the generalized Fried- 
mann equation. Indeed, when the electromagnetic terms and those measuring anisotropy are 
removed, the above assumes its familiar FRW form. Note that, on using (I4.3.7p . the Gauss- 
Co dacci formula (see Eq. (14.3.61) ) reads 

111 

K-ab = ~ TUi ah + E ab + - (7r ab + V a b) - - @(o" a fe + uj ab ) + a c{a a c b) - u c{a u c b) + 2cr c[a u c b ] . (4.3.8) 

The latter may be used to calculate the curvature of the 3-space along a chosen direction. For 
instance, TZ ab B a B b gives the 3-curvature distortions along the magnetic field lines. 



5 Non-linear cosmological magnetohydrodynamics 

With the exception of any period of inflation and early reheating, the universe has been a 
good conductor throughout its lifetime. As a result, B- fields of cosmological origin must have 
remained frozen into the expanding cosmic fluid during most of their evolution. This allows 
us to study the magnetic effects on structure formation within the ideal-MHD limits. 



5. 1 Ideal MHD approximation 

Consider a general spacetime filled with a single barotropic fluid of very high conductivity. 
Ohm's law (see Eq. (I3.4.2P ) guarantees that in the frame of the fundamental observer the 
electric field vanishes despite the presence of non-zero currents. Therefore, in the ideal MHD 
limit the energy-momentum tensor of the residual magnetic field simplifies to [99] 

T { ab ] = \ B 2 u a u b + \ B 2 h ab + U ab , (5.1.1) 
2 6 

with 

IU = n (o&> = ^ B 2 h ab - B a B b . (5.1.2) 

Accordingly, the 5-field corresponds to an imperfect fluid with energy density ps = B 2 /2, 
isotropic pressure pb = B 2 /6, and anisotropic stresses represented by the symmetric and 
trace- free tensor n a fe. 

Similarly, in the absence of an electric field, Maxwell's equations reduce to a single propagation 
formula, namely the covariant magnetic induction equation, 

B {a) = (a ab + e abc u c - 6/i afe ) B b , (5.1.3) 

and the following three constraints 
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car\B a = J a -e abc A b B e , (5.1.4) 
u a B a = ±p e , (5.1.5) 
D a S a = 0. (5.1.6) 

The right-hand side of (15.1.31) is due to the relative motion of the neighbouring observes and 
guarantees that the magnetic field lines always connect the same matter particles [99]. This 
means that the field remains frozen-in with the highly conducting fluid. Expression (15.1.41) 
provides a direct relation between the spatial currents, which are responsible for keeping 
the field lines frozen-in with the matter, and the magnetic field itself (e.g. see [H]). Note 
Eq. (15.1.51) which shows that rotating neighbouring observers will measure a non-zero charge 
density, triggered by their relative motion, unless uj a B a = 0. Finally, (15.1.61) demonstrates 
that in the absence of magnetic monopoles the field lines remain closed and B a is a solenoidal 
vector. 



5.2 Magnetic evolution 



The magnetic induction equation also provides the non-linear evolution law for the energy 
density of the field. More precisely, contracting Eq. (I5.1.3P with B a and then using (15.1.21) we 
arrive at 

(B*)' = ~QB 2 -2a ab W b . (5.2.1) 

This shows that in a highly conducting cosmic medium we have B 2 oc a -4 always unless there 
is substantial anisotropy, in which case the 5-field behaves as an anisotropic radiative fluid. 
In fact, in a homogeneous and anisotropic radiation-dominated universe this latter situation 
arises even when the anisotropy is small because close to isotropy U ab oc p and the evolution 
of B 2 is determined at second-order with B 2 /p ra d — > a/Q oc 1/ log(t) during the radiation era 
(see j6T|65] and also [66ll02j). 

The nonlinear evolution of the anisotropic magnetic stresses comes from the time derivative 
of (15.1.21) . which by means of Eqs. (15.1.31) and (15.2.11) leads to 

iU = -- en afe - ~ B 2 a ab + 2a c{a W b) - 2u c{a U c b) . (5.2.2) 



5.3 Conservation laws 



The energy momentum tensor corresponding to a magnetic single perfect fluid of infinite 
conductivity is given by 

T ab =(p+\ 5 2 ) u a u b +{p+l B 2 ^j h ab + U ab , (5.3.1) 
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Accordingly, the medium corresponds to an imperfect fluid with effective density equal to 
p + B 2 /2, isotropic pressure given by p + B 2 /6, zero heat flux and solely magnetic anisotropic 
stresses represented by U ab (see ( 15.1.21) ). 

When applied to the above, and using the MHD form of Maxwell's equations, the standard 
conservation law \7 b T ab = decomposes into the following expressions that respectively de- 
scribe the energy-density 

P = -(p + p)® (5.3.2) 

and the momentum-density 

[p + p + jj B 2 ^j A a = -D a p - e abc B b cni\ B c - U ab A b . (5.3.3) 

conservation [9911101] . Note the absence of magnetic terms in Eq. (15.3.21) . This is guaranteed by 
the magnetic induction equation (I5.1.3P and reflects the fact that the magnetic energy density 
is separately conserved|jj Also, when there are no pressure gradients, (15.3.31) gives A a B a = 
to ensure that the field exerts no forces along its own direction. Finally, the left-hand side of 
(15.3.31) shows that the magnetic contribution to the total inertial mass of the system is 2B 2 /3. 

Contracting Eq. (15.3.31) along B a we find that the contribution of the -B-field to the momentum 
density vanishes, thus guaranteeing that the magnetic Lorentz force is always normal to the 
field lines. Also, the second term in the right-hand side of (I5.3.3P decomposes as 

£ a6c 5 6 curl B c = i D a B 2 - B b D b B a . (5.3.4) 

The last term in the above is the result of the magnetic tension. In so far as this tension stress 
is not balanced by the pressure gradients, the field lines are out of equilibrium and there is a 
non-zero Lorentz force acting on the particles of the magnetic fluid. 

5.4 Magnetic tension 

The anisotropic nature of magnetic fields is encoded in the energy-momentum tensor of the 
field - and particularly in the anisotropic pressure tensor (see (I5.1.2P ). This ensures that the 
-B-fields exert a positive pressure orthogonal to their own direction, while carrying a tension 
along B a . Both of these very well known magnetic features are reflected in the eigenvalues of 
U ab , which are positive (1/3) perpendicular to B a and negative (—2/3) parallel to it. 

The magnetic tension also demonstrates the elasticity of the field lines and their tendency 
to remain as straight as possible by reacting to any effect that distorts them from equilib- 
rium [i~3][T6] . Within the context of general relativity, we can see this tendency by looking 

3 Expressions (|5.3.2|) and (|5.3.3p can also be obtained from the conservation laws (|3.3.ip and (|3,3.2p . 
This is done by taking the ideal-MHD limit of the latter, using Eq. (|5.1.4|) and assuming perfect-fluid 
matter. 
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at the effect of the field on the spatial curvature of a magnetic spacetime. In particular, 
contracting (14.3.61) along B a twice, and ignoring all sources but the magnetic field, we find 
that [H2] 

K ab B a B b = l -B* + ~ U ab B a B b = , (5.4.1) 

given that H a bB b = — (2B 2 /3)B a (see (15.1.21) ). Accordingly, despite the magnetic presence and 
the energy density contribution of the field, the curvature of the 3-space in the direction of the 
magnetic force lines is zero. Mathematically speaking, it is the contribution of the negative 
magnetic pressure along B a which cancels out the positive input of the field's energy density. 
More intuitively, however, one could argue that it is the elasticity of the magnetic lines, and 
their tendency to remain straight, that maintains the zero curvature along B a . 

In the presence of sources, 1Z ab B a B h is generally non-zero and the magnetic lines are forced out 
of equilibrium. When allowing for spatial inhomogeneities, the reaction of the field's tension, 
as expressed through the second stress in the right-hand side of (15.3.41) . to these geometrically 
induced distortions generally leads to counter-intuitive effects. The potential implications of 
these tension stresses, for cosmology as well as astrophysics, have been discussed in |143j - 
[146] . In spatially homogeneous spacetimes on the other hand, like the Bianchi I models, 
the magnetic Lorentz force vanishes and the field's tension is manifested only as a negative 
energy-density input, through the anisotropic pressure tensor (Il ab - see Eq. (15.1.21) ). Such 
cosmologies were recently investigated in [89"] . 



5.5 Kinematical evolution 



The magnetic presence affects the kinematics of the highly conducting medium both directly 
and indirectly. In particular, Raychaudhuri's equation reads 

® = -\® 2 -\(p + Zp + B 2 )- 2{a 2 - cu 2 ) + B a A a + A a A a , (5.5.1) 

showing that the field's contribution to the total gravitational mass is B 2 . There are also 
indirect magnetic effects propagating through the rest of the terms in the right-hand side of 
(15.5.11) . Of these effects, probably the most important are carried by the acceleration terms 
which follow from Eq. (15.3.31) . Similarly, there are direct and indirect magnetic effects on the 
shear evolution. In the magnetic environment the propagation equation of the shear tensor 
takes the form 

2 1 

o-{ab) = -g @a-ab - o- c ( a a c b ) - uj{ a uj b ) + D {a A b) + A {a A b) - E ab + - Il ab , (5.5.2) 

which reveals how the anisotropic pressure of the field directly induces and also affects kine- 
matical anisotropics. Rotation, on the other hand hand, is affected only indirectly through 
the magnetic effects on the 4-acceleration of the fluid. To be precise, 

2 1 

U( a ) = - 75 curl A a + a ab u . (5.5.3) 
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This means that the B-field not only affects the rotational behaviour of the observers' frame 
but it can also generate vorticity. 

The kinematical constraints are also only indirectly affected by the field's presence. This is 
manifested by the absence of explicit magnetic terms in Eqs. (j2.3.5j) - (l2.3.7p . When applied to 
a highly conducting barotropic fluid expressions (12.3.61) and (12.3.71) remain unchanged, while 
( 12.3.51) reduces to 

D 6 a afe = jj D a 9 + curl cu a + 2e abc A b u c . (5.5.4) 

When D a cx a b = 0, the shear is divergence-free and its pure-tensor part has been isolated. In 
other words, there is no way of using a a b to construct a scalar or a vector. This condition will be 
imposed later, in § [IT], when studying the magnetic effects of gravitational- wave perturbations. 



5.6 Spacetime curvature 



The energy-momentum tensor of a highly conducting magnetic perfect fluid is given by 
Eqs. ( 15TTTD . (15X21) . In that case formulae (HT5|I - (HT51> reduce to 



R ab u a u b = - (p + 3p + B 2 ) 



h a b R bc u c = 



(5.6.1) 



and 



h a c h b d R cd 



p-p+-B 2 



hab + n Q 



(5.6.2) 



respectively. We note that the above are fully non-linear expressions and that the former 
quantifies the total gravitational mass of the MHD fluid (in agreement with (15.5. ip ). 



We evaluate the magnetic contribution to the local 3-Ricci tensor orthogonal to the u a - 
congruence, by taking the MHD limit of the geometrical relations given in § 14.31 Expressions 
( 14.3.61) and (14.3.71) . in particular, simplify to 



2/1 1 

Tl ah = E ab + - [p + - B 2 - - 6 2 + a 2 - uj 2 ) h ab + £ U ab - ^ Q(a ab + uj, 



ab 



and 



3 V 2 3 J 2 3 

+c r c( a cr c b) - ^ c ( a uj c b) + 2a c [ a u c b] (5.6.3) 



V - 2{p+^B 2 -^Q 2 + a 2 -cu 2 ) . (5.G.1) 



respectively. This is simply the familiar Friedmann equation generalized to the magnetic and 
highly conducting environment. 
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5. 7 Evolution of inhomogeneities 



5. 7. 1 Basic variables 

Spatial inhomogeneities of a given physical quantity are covariantly described by its orthog- 
onally projected gradients. For our purposes, the basic variables are the comoving fractional 
gradients, 

A a = -D a p and B a = B a B 2 , (5.7.1) 

p B z 

describing spatial variations in the fluid and the magnetic energy density distributions (see [147] - 
|149j ) and [551100] respectively). Both variables are dimensionless by construction, and vanish 
identically when the spacetime is spatially homogeneous. As we will show later, A a and B a 
respectively describe spatial variations in the fluid and the magnetic energy densities, as mea- 
sured by a pair of neighbouring fundamental observers, in a gauge-invariant way (see § 16.3. ip . 

The above are supplemented by a pair of auxiliary variables, which describe spatial inhomo- 
geneities in the average expansion and the isotropic pressure of the fluid. These are [147] 

Z a = aD a Q and Y a = D a p, (5.7.2) 

respectively. It should be noted that, for barotropic fluids, the variable Y a is directly related to 
A a and therefore is redundant for all practical purposes. Finally, one may also monitor spatial 
inhomogeneities in the distribution of the magnetic field vector by means of the orthogonally 
projected tensor [99] 

M ab = B b B a . (5.7.3) 

This is trace-free in the ideal MHD limit (see constraint (I5.1.6P ) but, contrary to the gradients 
(15.7.11) and (I5.7.2p . it does not vanish unless the FRW metric is spatially flat (see Appendix A 
in [100] and also § 16.3.21 here) . 



5.7.2 Evolution equations 

In the presence of magnetic fields, the non-linear evolution of spatial inhomogeneities in the 
density distribution of a single, highly conducting perfect fluid is described by the expression 



A {a) =wQA a - (1 + w)Z a + — e abc B h cm\B c + - t£(l + w)aQA a 



a —e abc B b cni\B c +^ 
p 3 

- (o- ba + u ba ) A b + — Ii ab A b , (5.7.4) 



P 



with w = p/p and c\ = B 2 /p(l + w) representing the Alfven speed (see § 16.1.41 below). 
This equation has been obtained by taking the proper-time derivative of (15.7. la ) and then 
projecting it onto the observer's 3-dimensional instantaneous rest space. In the process, we 
have also employed the conservation laws (15.3.21) and (I5.3.3p . 
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Similarly, starting from definition (15. 7.2a ) and using (I5.3.3p . (15.5.11) . we arrive at the following 
non-linear evolution equation for the expansion gradients 



Z {a) = - 2 -®Z a - X -p^ a - X - B 2 B a + j| ae abc B b cni\B c + aB a A 



+ 



-n-3(cr 2 -uj 2 ) + A + A b A b 



aA a + 2aA b B a A b 



{a ba + uj ba )Z b + - aU ab A b - 2aD a (a 2 



(5.7.5) 



where A = D a A a and 1Z is the Ricci scalar of the observer's local 3-D rest space. 



Finally, the orthogonally-projected time derivative of ( 15.7.1b ) leads to the non-linear propa- 
gation formula monitoring spatial inhomogeneities in the magnetic energy density 



B 



(a) 



3(1 + w) 



4wG 
3(1 +w) 



4a6 



3p(l + w] 



£ nbr B b cui\ B c 



\aQ (l + lcl)A a 



- (a ba + uj ba ) B + 



3(1 + w) 



(a ba + u ba ) A b - 



4a9 



. U ab A b - -4 n bc D a (x bc 
3p(l + w) B 2 



~ a bc DJl bc + A a bc U bc B a - ^ a bc U bc A a . 



(5.7.6) 



In deriving the above, we have employed the non-linear relations (15.1.31) and (15.2.11) . as well as 
(15.7.41) . This helped to eliminate the expansion gradients from the right-hand sideof Eq. (I5.7.6p . 
We emphasize that (15.7.61) cannot beused when the equation of state of the cosmic medium 
has the vacuumenergy form contributed by a cosmological constant in the Einsteinequations, 
namely for w — p/p — — 1. In that case we can no longeruse (15. 7.4ft to eliminate Z a from (15.7.61) 
and the evolution of density inhomogeneities is described by a different system of equations 
(see § ELgp . 



6 Perturbed magnetic FRW cosmologies 



The symmetries of FRW spacetimes cannot naturally accommodate generic anisotropic sources 
like magnetic fields. This means that when based on a Friedmann background, any magnetic 
study will require some degree of approximation. Nevertheless, it is intuitively plausible that 
sufficiently weak _B-fields can be adequately studied within perturbed almost-FRW models. 
This belief, which has been at the core of almost every study of cosmological magnetic fields, 
has been confirmed by the analysis of exact [65] and perturbed magnetic Bianchi I uni- 
verses |102j . 
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6.1 Background evolution 



6.1.1 Zero-order approach 

We consider a spatially flat Friedmann background and allow for the presence of a weak 
magnetic field. This is effectively a test field that does not disrupt isotropy and has B 2 <C p. 
The background B-field can be either homogeneous or sufficiently random, so that (B a ) = 
7^ (B 2 ) (with (B 2 ) homogeneous on all the scales of interest). In both cases the magnetic 
anisotropic stress is treated as a first-order perturbation and the only zero-order magnetic 
variable is B 2 . First-order perturbations arise through H a b, the spatial gradients and the 
spatial derivatives DbB a - If the background test field is homogeneous, the scalar B 2 describes 
its local density and pressure. For a stochastic test field, B 2 describes the average density and 
pressure of the background field. In either case, the form of the linearized MHD equations 
guarantees that the magnetic effects propagate only via 5 2 -terms (see § 17.21) and that, at least 
when studying scalar (i.e. density) perturbations, the two approaches lead to the same linear 
equations and results. 

Alternatively, one can treat the magnetic field as a first-order perturbation, usually as a 
stochastic Gaussian field (e.g. see (76.80J). One can adapt our analysis to this scheme by 
removing the magneto-curvature terms from all the linear formulae (see § I7.2p . 



6.1.2 Zero-order equations 

Locally, the evolution of the background is governed by the standard Friedmann and Ray- 
chaudhuri equations. Introducing H = a/a as our zero-order Hubble parameter, these are 

3H 2 = p and H = -H 2 - - p(l + 3w) , (6.1.1) 

6 

and are supplemented by the conservation laws for the matter and the magnetic energy den- 
sities, namely by 

p=-3Hp(l + w) and (b 2 ) ' = -AHB 2 , (6.1.2) 

respectively^ This is obtained from the magnetic induction equation (see (15.1.31) ). which 
guarantees that 

B a = -2HB a . (6.1.3) 

Note that the weakness of the 5-field allows us to neglect the magnetic contribution to the total 
energy density and gravitational mass of the system. This explains the absence of magnetic 
terms from the right-hand sides of (16.1.1b ) and ( 16.1.1b ). provided w ^ —1/3 in the latter case. 



4 The spatially averaged magnetic energy density depends only on time and evolves in tune with its 
local counterpart. Indeed, using the local result (|6.1.2b ). it is straightforward to show that (B 2 )' = 
—AH(B 2 ) on any background domain because (H) = H in FRW models (e.g. see [150J for details). 
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Also, (16. 1.2b ) guarantees the radiation-like evolution B 2 oc a 4 for the zero-order field, which 
means that the background magnetic flux is conserved. 



6.1.3 Specifying the magnetic medium 

The nature of the cosmic medium is determined by its equation of state. Here, we will be 
dealing with a highly conducting barotropic fluid with p = wp, where w is the barotropic 
index. This vanishes in the case of dust, takes the value w = 1/3 for radiation and indicates a 
stiff fluid when w = 1. Also, for w < —1/3 the total gravitational mass of the matter is negative 
and the model enters a period of accelerated expansion, which corresponds to exponential de 
Sitter-type inflation when w = — 1. We find that 

w = -(1 + w)(c 2 s - w)e , (6.1.4) 

where c 2 s = dp/ dp is the square of the adiabatic sound speed. Therefore, a time-invariant 
barotropic index in a non-static universe (i.e. w = and 0^0) means that c 2 = w = constant, 
and vice versa. 



6.1.4 Magnetic strength and Alfven speed 



The magnetic strength, relative to the matter component, is measured by the dimensionless 
energy-density ratio (3 = B 2 / p. Then, given the weakness of the field (i.e. since j3 -C 1), the 
square of the Alfven speed is conveniently defined as 



P 



1 + w' 



(6.1.5) 



ensuring that C 1 at all times. The only exception is when w — > —1, in which case the 
proper definition for the Alfven speed is c 2 = Thus, from now on, w ^ —1 

unless stated otherwise. On using (I6.1.4P and Eqs. (16.1.21) . the time derivative of the above 
gives 

(c 2 J = -(l-3c 2 s )Hcl (6.1.6) 

to zero order. Accordingly, the Alfven speed does not change with time when c 2 = 1/3, which 
agrees with the radiation-like linear evolution of the magnetic energy density in almost-FRW 
environmentsOD 



5 In the standard literature the square of the Alfven speed is c 2 = B 2 / p = (3 (e.g. see |132| ). 

6 When shear anisotropy is included, the parameter and therefore c 2 are no longer constant during 
the radiation-dominated epoch but display a slow, 'quasi static', logarithmic decay (see [61|65j and 
also |66llU2j h 
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6.2 Linear evolution 



We linearize the full equations by treating quantities with non-zero background values as 
zero-order variables in perturbative terms. On the other hand, quantities that vanish in the 
background are of order one and higher-order terms are ignored. Therefore, the only zero-order 
quantities in our case are the fluid energy density (p) and pressure (p), the magnetic energy 
density (-B 2 ), and the Hubble parameter (H). 

The linear conservation law of the magnetic energy density maintains its background form, 
given by Eq. (16. 1.2b ). Similarly, the linear propagation of the magnetic vector is still given by 
(16.1.31) . which implies that B a oc a~ 2 and that the magnetic flux remains conserved. On the 
other hand, the scalar f? a curl B a , which is used to measure the helicity of the field (e.g. see [ZED, 
obeys the linear evolution law 

(B a cm\B a )' = -5HB a cm\B a . (6.2.1) 

The above means that the magnetic helicity decays as a" 5 and therefore it decreases faster 
than the energy density of the I?-field (recall that B 2 oc a -4 throughout the linear regime). 
Note that in deriving (16.2.11) we have assumed overall charge neutrality and used Eq. (I6.1.3p . 
together with the first-order expression 5 a (curl B a )' = —3HB a curl B a , obtained by means of 
commutator (1A.2j) in § IA.2I 



To first order, the conservation law of the fluid energy- density retains its the background form 
(see Eq. (16.1.2a )). However, when linearized around a FRW background, the momentum- 
density conservation equation (I5.3.3P reduces to 

(1 + w) ( 1 + ~ cf) pA a = -c 2 B a p - e abc B b cm\ B c , (6.2.2) 



3 

with D a p = c 2 D a p due to the barotropic nature of the medium. 

The magnetic field's presence affects the linear kinematics in a number of ways. The average 
volume expansion, for example, is determined by the following version of the Raychaudhuri 
equation 

11 r 2 / 2 \ r 2 1 



O = — O 2 - - p(l + 3u; + 0) - , s r 1 - - < D a A a - D a B a + -ciTZ, (6.2.3) 
3 2 HK H! a{l+w)\ 3 V 2a 3 a V ; 

where 1Z = h ab 7Z a b is the 3-Ricci scalar and we have kept up to c 2 -order terms. As expected, any 
local increase in the magnetic energy density will decelerate the linear expansion of the model 
in the same way that matter perturbations do. However, the 5-field has an additional effect 
due to the magneto-curvature term in the right-hand side of (I6.2.3p . This further decelerates 
the expansion if 1Z < 0, but tends to accelerate it when 1Z is positive. This counter-intuitive 
behaviour, whichresults from the tension of the field, can have nontrivial and veryunexpected 
implications for spatially curved magnetic spacetimes |143j -[146j. 
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The rotational behaviour of the weakly magnetic fluid is described by the vorticity propagation 
equation, which to linear order takes the form 



-2 



2 s 



1 



c?|- 



6 



1 



HUn 

2p(l 



W] 



1 



" 2X B h D h cm\B n 



(6.2.4) 



This means that the magnetic presence will generate vorticity unless the field's distribution 
is curl-free, or more specifically, unless curl B a remains unchanged along the direction of the 
magnetic force lines. Thus, I?-fields can act as sources of rotation. The magnetic effect on 
pre-existing vorticesis more difficult to quantify. As we will see later, in section§ 110.31 the 
presence of the field can slow down the standard decay rate of vorticity and therefore increase 
the residualamount of rotation relative to magnetic-free cosmologies. Recallthat in the absence 
of the B- field, vorticity always decays unlessc^ > 2/3 [151J. 



The evolution of linear kinematical anisotropies in the magnetic fluid is governed by the 
following expression of the shear propagation equation 



1 „ „ ct (. 2 2 \ „ . cl 



a ab = -2Ha ah + - U ab - E ah - a ^ s +w) [1 ~ g <) D{ a A 6> - ^ D {a B b) 

+^^y (l - \ci) B c D c D {a B b) + ±<*K {ab) . (6.2.5) 

The magnetic effects are diverse. The anisotropic pressure of the field is supplemented by 
anisotropies in the distribution of the magnetic vector and by those in the energy density of 
the B- field. In addition, there is a purely geometrical magneto-curvature effect, provided that 
Hub) i s non-zero at the linear level. 

In the ideal MHD limit the kinematical variables also satisfy constraints (15.5.41) . (I2.3.6P and 
(12.3.71) . which then linearize to 

D 6 (7 o6 = ^D a e + curlu; a , B a u a = (6.2.6) 

and 

H ab = curl a ab + B {a u b) , (6.2.7) 

respectively. According to the above, uj a is a linear solenoidal vector and in the absence of 
rotation the magnetic Weyl component is fully determined by the shear tensor. Additional 
(geometrical) constraints between the kinematical and the dynamical quantities in the highly 
conducting medium are provided by (I5.6.3|) and (15.6.41) . The respective linearized counterparts 
of these expressions are 

n ab = E ab + - (p+\ b2 -\ 02 ) h ab + 2 n ^ - H (°ab + "ah) (6.2.8) 

and 

TZ = 2[p+ l -B 2 - l -Q^ . (6.2.9) 
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6.3 Spatial inhomogeneities 



6. 3. 1 Local interpretation of the inhomogeneity variables 

To linear order, the inhomogeneity variables defined in section 15.7.11 describe measurable 
differences in the spatial distribution of physical quantities. In order to verify this, consider the 
relative position vector Xa connecting the same two points on the worldlines of neighbouring 
fundamental observers. Following [12211147] . we have XaU a = and 

Xa = Hxa + (cr ab + u ab )x h , (6.3. 1) 

which means that in a FRW universe Xa — a>X° a (with x° a = constant and = 1). If P and P 
are points on the two worldlines with coordinates (x a ) and (x a + Xa) respectively and p, p are 
the associated values of the energy density, a Taylor expansion around P gives 

p-p = X a V a p, (6.3.2) 

to leading order. Since Xa oc a in the FRW background and using definition (15. 7. la ), the above 
translates into 

Sp = x a A a , (6.3.3) 

where bp — (p — p)j p. In other words, the comoving fractional gradient A a describes the mea- 
surable local density variation between two neighbouring fundamental observers. Moreover, 
A a closely corresponds to the familiar energy-density contrast of the non-covariant studies. 
The same analysis also applies to the rest of the variables defined in 15.7.11 



6.3.2 Gauge invariance of the inhomogeneity variables 

Cosmological perturbations have long been known to suffer from gauge-related ambiguities, 
reflecting the fact that in perturbation theory one deals with two spacetime manifolds. The 
first is the physical spacetime (W) that corresponds to the real universe, while the second 
(W) is a fictitious background described by an idealised mathematical model. To proceed 
we need to establish an one-to-one correspondence, namely a gauge (f> : W — » W, between 
these two manifolds. Such point identification maps are generally arbitrary and the gauge 
problem stems from our inherent freedom to make gauge transformations. The latter differ 
from ordinary coordinate transformations because they also change the point identification 
between W and W. 

We define perturbations as the difference between the value of a given quantity at some event 
in the realistic universe and its value at the corresponding (through the gauge) event in the 
background spacetime. Gauge transformations, however, generally change this correspondence 
and therefore the value of the perturbation itself. This makes perturbations gauge- dependent 
and arbitrary. For instance, we can select the gauge in such a way that a given perturbation 
vanishes [147] . Following the Stewart & Walker lema, the simplest quantities that remain 
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invariant under gauge-transformations are scalars that are constant in the background space- 
time and tensors that vanish there [152] Note that the same general criteria also apply to 
second-order perturbations, but this time the Stewart & Walker requirements must be satisfied 
by the first-order variables [153J. 

In a Friedmann universe all physical variables are functions of cosmic time only, though they 
acquire additional spatial dependence in perturbed 'nearly-FRW universes. As a result of 
the spatial homogeneity of the unperturbed Friedmann background, the variables (15.7.11) and 
( 15.7.21) vanish to zero order. Consider the density gradient for example. Then, 

D a p = h a b V bP = h a °V p + h a Q V aP = , (6.3.4) 

because V a p = and in a comoving frame h a ° = (with a = 1, 2, 3 - see also P-47J). This result 
guarantees the vanishing of the inhomogeneity variables (15.7.11) . (15.7.21) . in the background, and 
therefore ensures their gauge invariance at the linear perturbative level [152J. Note, however, 
that the gauge invariance of D b B a , and therefore of A4 ab , are guaranteed only in spatially flat 
FRW backgrounds [100] . One can see that by applying the 3-Ricci identity to a zero-order 
3- vector v a (see Eq. (IA.2[) in § lA.ip . In the absence of rotation the latter reads 

2D [a D b] v c = 1Z abcd v d , (6.3.5) 

which means that D b v a ^ unless lZ abc d = [100J. 



6.3.3 Irreducible inhomogeneity variables 

The orthogonally projected density gradient, A a , contains collective information about three 
types of inhomogeneities, namely density perturbations, vortices and shape-distortions. This 
information is encoded in the dimensionless comoving gradient 

A q6 = aD b A a (6.3.6) 

and can be decoded by splitting the latter into its irreducible components as follows [149|I154] 

A ab = - aD c A c h ab + aD [b A a] + aB {b A a) . (6.3.7) 

The quantities on the right-hand side are associated with the aforementioned three different 
types of inhomogeneity. The scalar D a A a describes spatial variations in the matter density 
(i.e. overdensities or underdensities), the skew part is related to magnitude-preserving changes 
of A a (i.e. rotations), and the symmetric and trace-free tensor D/ b A a \ describes shape distor- 
tions in the anisotropy pattern of the gradient field (e.g. pancakes or cigar-like structures). 



7 The only alternative way of constructing gauge-independent variables is by using tensors which 
are linear combinations of products of the Kronecker deldas with constant coefficients |152] . 
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In an exactly analogous way, all the information regarding perturbations in the volume ex- 
pansion and in the magnetic energy density is stored in the dimensionless second-rank tensors 

Z ab = aD b Z a and B ab = aD h B a , (6.3.8) 

respectively. These decompose into their irreducible parts according to [149 j 101] 

Z ab = ^ aD c Z c h ab + aD [b Z a] + aD {b Z a) (6.3.9) 

and 

B ab = ^ aD c B c h ab + aD [b B a] + aB {b B a) . (6.3.10) 



6.3.4 Linear equations 

When linearized around a FRW background permeated by a weak magnetic field, the non- 
linear propagation equations of section 15.7.21 reduce to 



3aH 

A a = 3wHA a - (1 + w)Z a + e abc B b cvn\B c + 2c£(l + w)aHA a , (6.3.11) 

P 

Z a = -2HZ a - ~ pA a - ~ B 2 B a + ~ ae abc B b cm\ B c + aD a A (6.3.12) 

and 

= 57TZ S A « - 7T~~" A " " n ^ \ ^ S ' curl ^ C " 4a # f 1 + I c ') ^ , (6.3.13) 
3(1 + w) 1 + w p(l + w) V 3 / 

respectively. Therefore, the linear evolution of the inhomogeneities depends on the 4-acceleration 
of the highly conducting matter. To first order, the latter satisfies the momentum-density con- 
servation law (see Eq. (I6.2.2P ) 

(1 + w) (l + ^ tity a P A a = -c 2 sP A a - ae abc B b cur\ B c . (6.3.14) 

Using the commutation law D[ a D b ]t> c = —u ab V{ c ) + T^dcbaV d /2 between the spatial gradients 
of v a (with v a u a = - see Eq. (1A.2I) in § IA.1I) ). and the constraint (I5.1.6p . the projected 
divergence of the above leads to 

where 

n = 2[p+ l -B 2 - l -Q^ , (6.3.16) 
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is the linearized 3-Ricci scalar (see Eq. (I6.2.9P ). Substituting (16.3. 14[) into the right-hand side 
of (16.3.131) we arrive at 

4 • 4(c 2 -w)H 
B a = A a + 4lC ; W)ti A a , (6.3.17) 

3(1 + w) 1 + w 

which means that linear inhomogeneities in the magnetic energy density are only affected 
by those of the matter. To close the system of (16. 3. lip . (16.3.121) and (16.3.171) we require the 
propagation equation of the 3-Ricci scalar. Taking the time derivative of (I6.2.9P and linearising 
we obtain 

1l = -2 (l + \ cl) HK + n 4 f A HD a A a + ^ B a B a . (6.3.18) 

V 3 / ail + w) V 6 / a 



We note that for c 2 = w, as it happens during the radiation and the dust epochs when the 
equation of state of the cosmic medium remains unchanged (i.e. when w = - see Eq. (I6.1.4P ). 
the above expression reduces to 

A - ( 6 - 3 - 19 ) 

3(1 + w) 

This guarantees that linear perturbations in the magnetic energy density evolve in step with 
those in the density of the highly conducting matter. When dealing with a radiative fluid the 
above gives B a = A a , while B a = 4A a /3 for non-relativistic dust{f] The implication of (16.3.191) 
is that the effective entropy perturbations coming from the different dynamical behaviour of 
the two components (i.e. the fluid and the magnetic field) are constant, or zero (see section [9]). 



7 Density perturbations 



7.1 Basic variables 



Decomposition (I6.3.7P allows for a fully covariant and gauge-invariant study of linear density 
perturbations, vortices and shape distortions. Here, we will focus on the scalar component of 
A a b. This is defined by 

A = aD a A a , (7.1.1) 

and when positive (negative) it describes overdensities (underdensities) in the distribution 
of the perturbed fluid. In other words, A corresponds to the relative increase or decrease in 
the matter density as measured by a pair of neighbouring fundamental observers. Similarly, 
Z = aD a Z a and B = aD a B a describe scalar perturbations in the volume expansion and the 
magnetic energy density respectively. 

8 The covariantly derived result (|6.3.19p is in complete agreement with relation B oc p 2 /^ l+w \ 
familiar from conventional studies of weakly magnetized FRW models. Indeed, taking the square of 
the above, and then its variation, we find that 5B 2 /B 2 = A5p/3p(l + w) with 5B 2 /B 2 and Sp/p 
corresponding to B a and A a respectively. 
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To linear order definition (I7.1.ip reduces to A = (a 2 /p)D 2 p, while analogous expressions hold 
for B and Z. In addition, it helps to introduce the following rescaling 

JC = a 2 1l, (7.1.2) 

of the perturbed 3-Ricci scalar 1Z. This is also gauge invariant, given the spatial flatness of 
the FRW background. 



7.2 Linear equations 



Taking the orthogonally projected divergence of Eq. (16.3. lip , using (16.3. 14[) and keeping up to 
c 2 -order terms, given the weakness of the magnetic field, we arrive at the following equation 

A = 3w (l - ~ cf) HA — (1 + w)Z + | c 2 (l + w)HB - c 2 (l + w)HK , (7.2.1) 

for the linear evolution of A. Note that in deriving the above we have also set w = to 
zero order, which means that w = c 2 s = constant in the background, and employed the linear 
relation aD a A a = A. Following (17.2.11) . the field will generally act as a source for density 
perturbations even when there are no such distortions present initially. Also, the magnetic 
field's presence has a direct and an indirect effect on A. The former results from the pressure 
part of the Lorentz force (see decomposition (15.3.41) in § 15.31) and carries the effects of the 
isotropic magnetic pressure. The latter comes from the tension component of the Lorentz 
force and it is triggered by the magnetic coupling to the spatial curvature of the perturbed 
model. Surprisingly, a positive 3-curvature perturbation causes A to decrease, while a negative 
K, has the opposite effect. This rather counter- intuitive behaviour of the magneto-curvature 
term in ( 17.2. ip is a direct consequence of the elasticity of the field lines (see also Eq. ( 17.2. 2ft 
below). 

Similarly, using (I6.3.15p . together with the linear results aD a Z a = Z and aD a B a = B, the 
linearized orthogonally-projected divergences of ( 16.3. 12[) and ( 16.3. 17(1 lead to 



Z = -2 + HZ-\p(l-^A + \cl{l + w)pB- l -cl{l + w)pK 



c 2 



4 D'A - - c 2 J) 2 B (7.2.2) 



3 V 2 
and 

4 • 4(c 2 - w)H 
B= , r A+ [s A , (7.2.3) 

3(1 + w) 1 + w K ' 

respectively. Finally, starting from the linear propagation equation of the 3-Ricci scalar we 
obtain 

K = ~ <*HK + -^2- (l - % c 2 ) HA + 2<%HB . (7.2.4) 
3 a 1 + w \ 3 V a v 1 
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The system (I7.2.ip - (l7.2.4p describes the linear evolution of scalar inhomogeneities in the den- 
sity distribution of the matter in a weakly magnetic and spatially flat almost-FRW universe I" 9 ! 
When the cosmic medium has a time-independent equation of state (i.e. for w = 0) we have 
c 2 = w = constant and Eq. (17.2.31) reduces to 



B = - -A. (7.2.5) 

3(1 + w) V ; 



It should be noted that in [9"§]- |101| factors of the form lie 2 in the coefficients of the perturbed 
variables were set to unity, because of the overall weakness of the 5-field. Here, retaining the 
factors lie 2 in the system of (I7.2.ip - (I7.2.4I) has improved the accuracy of our linear equations 
and will help us refine the results of [99]- [101j . When 1 ± c\ ~ 1, expressions (I7.2.ip - (I7.2.4I) 
immediately reduce to their corresponding formulae of [99j- |101j . The reduction is obvious and 
straightforward for Eqs. (I7.2.1l) - (l7.2.3p but not for (17.2. 4p . There, one needs to use definition 
(17.1.21) to show that /C + (Ac 2 J?>)HfC ~ /C when 1 ± c 2 ~ 1. This requirement was overlooked 
in [103] , resulting in the inconsistent linearisation of their equations. 



7.3 Evolution in the radiation era 



During the radiation epoch the background dynamics is determined by the parameters w = 
1/3 = Cg, H — l/2t and p = 3/4t 2 . In addition, given the weakness of the magnetic field 
we have c 2 = 3/3/4, where f3 = B 2 /p = constant < 1. To proceed further we harmonically 
decompose the perturbed variables. In particular, we set A = J2k ^ k ^Q(k), where Q^) are the 
standard scalar harmonics (with k representing the associated wavenumber) and D a A( fc ) = 0. 
The harmonics are time-independent functions (i.e. Q^) = 0) and satisfy the scalar Laplace- 
Beltrami equation 




D 2 Q (k) = - - Q (k) . (7.3.1) 



Similarly, the rest of the perturbed variables decompose as Z = Z^Q^, B = B^Q^) and 
/C = IC^Q^k), with D a Z^ = = D a B^ = D a K,^ k \ Implementing this harmonic decomposi- 
tion, Eqs. dE2HD-dE£U now read 



9 When B 2 is treated as a first order perturbation, the 3-Ricci terms in the right-hand side of 
Eqs. (|7,2.ip . (|7.2.2p are second order and therefore vanish at the linear level. In that case linearized 
magnetized matter fluctuations evolve free of spatial curvature effects and (I7.2.4P is no longer re- 
quired. Note that c 2 B = a 2 D 2 B 2 /p by definition and therefore the variable c 2 B is still linear. 
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A (fc) = ±(i-lp) r 1 a« - - z<® - - pr x K.v> + ~ ^r^w , (7.3.2) 

2 \ 2 / 3 2 4 

2« = - f 1 + - ^ r 1 ^ - - (1 - /3)r 2 A<*> - - prW* + - (3t- 2 B^ 
V 2 / 8 8 16 

4 f ^ 2 f 1 _ 1 ^ A« + U-X (3B^ , (7.3.3) 



4 \a \ 2 J 8 \a 



t<® = ~ pr l K^ + i (1 - X - p) r J AW + ^ /sr 1 ^ , (7.3.4) 

£«=A (/c) . (7.3.5) 



During the radiation era the Hubble radius, in comoving proper time t, is given by A# = 
l/H = It. Also, physical wavelengths and comoving wave numbers are related by A/t = a/k. 
Substituting these expressions into Eq. (17.3.31) . and keeping up to /9-order terms (recall that 
P < 1), the system fl7X2l - fl7T33|) becomes 



AW = I ( i - 1 p) r 1 a« - - - - pr x KF> + ? /ft-^w 



i + -p)r 1 zW-~(i-p)r 2 



i / a 



8 lb 



(7.3.6) 



(7.3.7) 

(7.3.8) 
(7.3.9) 



where Xn/^k = k/aH. The above describe the evolution of linear matter perturbations in a 
weakly magnetic almost-FRW universe filled with a single highly conducting perfect fluid. 



7.3.1 Super-horizon scales 

When dealing with perturbations on super-Hubble lengths we have ^> A#. Then, given 
also that P -C 1, Eqs. fl7.3.6p - fl7.3.9l) reduce to the scale- independent system 

A = l{ 1 ~l P ) rl A " \ Z ~ \ ^ llC + \ Pt ~ lB ' (7 - 3 - 10) 

Z = -(l + \p) r x Z - | (1 - P)r 2 A - | /3t- 2 /C + A /3t- 2 i5 , (7.3.11) 

£ = -~ /ft" 1 * + I (l - I r 1 A + ^ /Sr 1 *? , (7.3.12) 

B = A. (7.3.13) 
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These can be solved analytically and lead to the following simple power-law evolution for the 
density contrast 

A = A(t) =C + dt z \ (7.3.14) 

where Co, Ci are constants (with i — 1, 2, 3) and the Zi are roots of a cubic algebraic equation. 
Ignoring terms of order 0((3 2 ) and higher, the latter reads 

8^ 3 - 4(1 - (3)z 2 - 4(1 - 0)z = 4(3 . (7.3.15) 

At the weak-field limit the above is solved perturbatively giving Z\ = — 1/2 + 5/3/6, z<i = —(3 
and z 3 = 1 — (3/3, which are all physically allowed solutions. Consequently, at this level of 
approximation, the linear evolution of A proceeds as a sum of power-laws: 

A = C + dt-\ + P + C 2 t~ p + C 3 t 1- ^ . (7.3.16) 

Note that in the absence of the magnetic field (i.e. for (3 = 0), we recover the standard 
evolution for the density contrast, familiar from the linear study of perturbed magnetic-free 
FRW models (e.g. see [155]- [i~57] ). Thus, in the weak- field limit the main magnetic effect is 
to reduce the growth rate of the dominant density mode. In addition, the field also decreases 
the rate of the standard decay mode and introduces a new 'non-adiabatic' decay mode] 10 



7.3.2 Sub-horizon scales 

Well below the horizon scale, A#/A& ^> 1 and the scale-dependent terms inside the brackets 
of Eq. (I7.3.7P become important. Then, on sub- horizon scales the system (I7.3.6l) - (l7.3.9p reads 



A« = -fl--(3) r x A^ - - - - (3r l K^ + - (3t~ l B^ , (7.3.17) 
2 \ 2 / 3 2 4 

-~ /3r 2 /c« + L p {**\ 2 t-ir^w , (7.3.18) 

W = -~ pr 1 ^ + i fi - I p) r x AW + j pr x BW , (7.3.19) 

B {k) =A^, (7.3.20) 

given that during the radiation era \n/^k = {Xh / \k)o{t/to) 1 ^ 2 ■ 



10 The magnetic suppression of the growing A-mode was first observed in [100] , although there the 
3-curvature effects were switched off. Allowing for the curvature effects, the analysis of |1U1| showed 
no magnetic effect on the growing mode to lowest order in (3, but a slight magnetic enhancement 
at higher order (see solution (54), (55) there). Here, the refined solutions showed a suppression 
proportional to (3 in line with the results of [100J. 
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When the 3-curvature effects are switched off, as is physically appropriate during the early 
universe, the system can be solved analytically. In particular, the solution for the linear matter 
perturbations reads 

A (») = A (*)( t ) =Cl+ + C 3 t 1/4 Y 1/2+ Mx) , (7.3.21) 

where are the integration constants (with i — 1, 2, 3) and Ji/2+/3(x) an d ^1/2+pix) are 
Bessel functions of the first and second kind respectively, with arguments determined by 

and c s = 1/VS. Also, recalling that /? <C 1 and keeping up to /5-order terms, we arrive at the 
following large-scale solution for B: 

B (k) = B (k ){t) = _2 ^ ^ _ 1 ^ rl + C2t i/4j i/2+/3(x) + c 3t i/4 Yl/2+/3 ( x ) . (7.3.23) 

As expected (see Eq. ( I7.3.20p ). the above is identical to (I7.3.2ip up to a constant. Setting 
A (fc) = A W and g(fc) = g(A) initia n y; so l u tions (17.3.211) and (I7.3.23D combine to give d ~ 
3/3(Af, fc) -^ fe) )/2. Then, 

A (fc) = 2/3 (A? ) - + C 2 t 1/4 J 1/2+ ,(x) + Cat 1 /^!/^^^) . (7.3.24) 

Based on the weakness of the magnetic field, we may approximate the two Bessel functions 
in the above given solutions by J1/2+/3 — J1/2 and Y1/2+P — ^1/2 respectively. In this case 
( I7.3.24P reduces to 




(7.3.25) 



where a = c s (Xh / ^k)oto (1 + /3/2). Therefore, we find that small-scale matter perturbations 
oscillate like magneto-sonic waves. The magnetic presence tends to reduce the amplitude of 
the oscillation and increase its frequency (in agreement with [lOOflOl] ). In both cases the effect 
of the field is proportional to its relative strength (i.e. to the ratio 3 = B 2 /p). As pointed out 
in [63], the increased frequency of A'*' should bring the peaks of short-wavelength oscillations 
in the density of the radiation component closer. This in turn could produce a potentially 
observable signature in the CMB. 

An additional magnetic effect arises from the presence of a constant mode in solution (17.3. 25p . 
This suggests that, unlike the magnetic-free case (e.g. see [156] ). magnetic matter perturba- 
tions in the pre-equality universe oscillate around a generally non-zero average value. This 
depends on the relation between A and B initially. For example, isocurvature initial condi- 
tions typically correspond to zero perturbation in the total energy density [16 5lll67j . In our 
case p t = p + B 2 /2 and the isocurvature requirement translates into A = (3/2)Bq (recall 
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that (3 = constant throughout this epoch). Following (17.3. 25|) . the latter implies that by the 
end of the radiation era the associated density contrast oscillates around the average value 
(A) ~ 3A . We finally note that, strictly speaking, the analysis of § 17.31 holds as long as 
the electrons are still relativistic (i.e. up to T ~ T e ). Subsequently, and until the end of the 
radiation era, one should examine the magnetic effects on perturbed non-relativistic matter 
on an expanding radiative background. 



7.4 Evolution in the dust era 



After the end of the radiation era, the unperturbed background is well approximated by 
w = = cj., H = 2/3t, p = 4/3t 2 . Also, c 2 A = (3 and it is no longer constant but decreases 
with time according to (3 oc a" 1 oc t~ 2 ^ 3 . Applying the usual harmonic decomposition to the 
perturbation variables, and keeping up to /3-order terms, the system (I7.2.1l) - (l7.2.4p reads 



3 



4 



4 



. 1 

3 

f -2 B (k) 



4 



A (fc) , 



(7.4.1) 

(7.4.2) 

(7.4.3) 
(7.4.4) 



with A# = 3t/2. The above system has no simple analytical solution because during the dust 
era the dimensionless parameter (3 is no longer constant. Nevertheless, as we shall see next, 
one can still extract useful qualitative and quantitative information about the magnetic effects 
on density perturbations in the post-recombination universe. 



7.4-1 Magnetically-induced Jeans' length 

The wave equation of the density contrast reveals the role of the magnetic pressure against 
gravitational collapse and helps to establish the scale of the magnetic Jeans' length. In par- 
ticular, taking the time derivative of (17.4.11) and using the rest of the formulae we arrive 
at 




t -2 A (k) 



+-(3r 2 K^. 
9 



(7.4.5) 
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Note that, in deriving the above, we have set B^ = 4A^ /3. This guarantees that B^ = 
4A( fc V 3 always (see Eq. (jTXl"]) ). 

The first term on the right-hand side of f I T . 4 . 5 1) shows that the magnetic presence reduces 
the diluting effect of the expansion by a small amount proportional to the field's relative 
strength. The second term gauges the opposing effects of gravity and magnetic pressure. 
When the curvature term is removed form (17.4.51) . and the quantity in the angled brackets is 
positive, gravity prevails and the density contrast grows. When the brackets take a negative 
value, the pressure of the magnetic field dominates and forces the inhomogeneity into an 
oscillatory phase. In the latter case the magnetic effect is relatively weak and also decays in 
time. Gravitational attraction and magnetic pressure balance each other out at the threshold 
1 — (8/3/3) [1 + (A#/Afc) 2 /6] = 0. This also determines an associated Jeans length, defined 
as the scale below which the magnetic pressure gradients dominate and the inhomogeneities 
oscillate. This 'magnetic' Jeans length is given by 



' ' Xh (7.4.6) 



J V 9 - 24/3 

and it is considerably smaller than the corresponding Hubble radius since (3 <C 1 (see also 
for analogous expressions). The anisotropy of the CMB constrains the current magnitude 
of a (homogeneous) cosmological magnetic field below 1CT 9 G [63f66] . In that case B$ ~ 
10 -58 GeV 4 , (3q ~ 10 -11 (recall that p = p cr ~ 10~ 47 GeV 4 today) and the associated magnetic 
Jeans length is 



Xj ~ KT 5 A H ~ 10 Kpc, (7.4.7) 

today. Random 5-fields, however, can have magnitudes up to 10~ 6 G in today's values. If 
B ~ 10~ 7 G, a strength supported by observations in galaxy clusters and high-redshift 
protogalactic structures, we find that 

Xj ~ HT 3 A H ~ 1 Mpc, (7.4.8) 

at present. This number, which can also be obtained through a fully Newtonian treatment |95j . 
determines the dimensions of the region that is gravitationally supported by the magnetic 
pressure gradients, and lies intriguingly close to the size of a galaxy cluster. 



7-4-2 Late-time evolution 

One can obtain an analytical solution for the evolution of magnetic density perturbations 
in the dust era by taking the time derivative of (I7.4.5p . Then, recalling that Xh/Xj. = 
(Xh/ Xk)o(t/toY^ 3 , using (17.4.31) and ( 17.4.41) . and keeping up to /3-order terms we arrive at 
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(7.4.9) 



The constant Cg is inversely proportional to the scale in question and vanishes when Bq = 
4A /3. As we will see in § the aforementioned initial conditions guarantee effective adia- 
baticity. Ignoring the weak and time-decaying terms from the right-hand side of (17.4.91) . we 
obtain the following late-time solution (see also Eq. (67) in [10 lj) 



A (fc) = Q^ai + ^2 + C 3t -2/3 + Q 



4 • 



with 
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(7.4.10) 



(7.4.11) 



In the absence of the 5-field (i.e. when (3q = 0), we immediately recover the standard magnetic- 
free solution with a% = 2/3 and «2 = — 1 (e.g. see [155] - [T57] ). Therefore, the main magnetic 
effect is to reduce the growth rate of density perturbations by an amount proportional to its 
relative strength (i.e. to the ratio /3 = (B 2 /p) ). It should be noted that the inhibiting role 
of the field was first observed in the Newtonian treatment of [90] arid later in the relativis- 
ts studies of |100|I101| . According to (17.4. 10p and (17.4. lip , the magnetic impact is inversely 
proportional to the scale in question. On super-horizon lengths, in particular, the above given 
solution reads 

A (fc) = Cit 2/3 + C2f -i + C3f -2/3 ^ (7.4.12) 

since C4 ~ (/3o/3)(\h/ \k)o[Bo — (4/3)A ] 1. Hence, on large scales the introduction of the 
.B-field simply adds the decaying t~ 2//3 mode to the standard magnetic-free result. Note also 
that A describes the directionally-averaged gravitational clumping of the matter. Generally, 
the perturbations will grow at different rates parallel and perpendicular to the magnetic field 
and so there will also be non-spherical evolution in the shapes of these distortions. 



7. 5 Evolution in false-vacuum environments 



False vacuum cosmological environments have been largely associated with inflation, since 
they typically lead to exponential, de Sitter-type expansion. Dynamically, such environments 
can be achieved by introducing a fluid with a non- conventional equation of state. Here, we will 
consider linear inhomogeneities in the energy density of a barotropic and of a non-barotropic 
magnetized medium with a false-vacuum equation of state. 

Nonlinear spatial inhomogeneities in the density distribution of a perfect, though not neces- 
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sarily barotropic, fluid propagate according tcQjjj 

A {a) = #A a - — D a p - (p + p)Z a - 3a(p + p)#A a - (a ba + u ba )A b , (7.5.1) 
P P 

The 4-acceleration vector seen above satisfies the momentum-density conservation law. In the 
absence of a magnetic field the latter reads 

(p + p)A a = -B a p, (7.5.2) 

with p+p representing the inertial mass of the fluid. When dealing with a magnetized medium, 
however, fl 7 . 5 . 2 1) is replaced by 



p + V + 3 B 2 j A a = -D a p - e abc B b cm\ B c - U ab A b , (7.5.3) 

which incorporates the magnetic contribution to the total inertial mass of the system. 

One way of achieving the false vacuum environment is by introducing an effective (barotropic) 
perfect fluid with p = —p. In that case Eq. (17.5.21) acts as a constraint demanding that 
DaP = 0. This immediately implies D a p = and therefore zero density inhomogeneities. This 
strict constraint is a direct and unavoidable consequence of having a medium with zero inertial 
mass. When a magnetic field is present, however, the 4-acceleration is given by (17.5.31) and 
this expression allows for nonzero pressure gradients even if p + p = 0. In that case we have 
D a p = — D a p ^ and Eq. (17.5.11) reduces to 

A a = 0, (7.5.4) 

at the linear perturbative level. Accordingly, all three types of linear density inhomogeneities 
(i.e. matter perturbations, vortices and shape distortions) remain constant as long as p+p = 0. 
Note that we always assume that B 2 <C p, which ensures that the magnetic field has no effect 
of the exponential expansion of the model (i.e. p + 3p + B 2 ~ p + 3p = —2p < 0). 

Standard de Sitter inflation is typically achieved by means of a single minimally coupled scalar 
field (tp), which generally does not behave like a barotropic medium. Instead, the <p-field acts 
as an effective perfect fluid with energy density and pressure given by 

p = p {Lp) = i <p 2 + V(ip) and p = p ( ^ = - (p 2 - V(ip) , (7.5.5) 

respectively (e.g. see |158|I159| ). Therefore, the scalar field corresponds to a barotropic fluid 



only when its kinetic or potential energies vanishes! 12 1 In the first instance we have (p = 
and p (v) = -p^) = -V{tp>), while in the second Eqs. ( 17.5. 5p give p^ = p^) = (p 2 /2. In any 
other case p = p — 2V((p) = (p 2 — p. The perfect-fluid description of a minimally coupled scalar 



11 One can derive (|7.5.ip from first principle, namely by taking the proper-time derivative of definition 
(|5.7.1r ). or directly from (|5.7.4p by substituting the 4-acceleration from Eq. (|3.3.2p . 

12 Non-minimally coupled scalar fields generally correspond to imperfect fluids with more complicated 
expressions for their associated effective density and pressure (e.g. see |160| ). 
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field refers to a family of observers moving with the (timelike) 4-velocity u a = —V a l fi/ ! f, 
where V a (p 7^ and <p = u a V a ip. Relative to the same observers, D a ip = and therefore 
D a p^ = D a p {Lp) = (pD a <p H5S]. This last result indicates a major change relative to the 
barotropic-fluid case and it will have a profound effect on the evolution of the inhomogeneities. 

The de Sitter phase occurs when the scalar field rolls slowly down its potential. Throughout 
this slow-rolling regime, (f 2 -C V(<p) and expressions (I7.5.5P imply that p^ 1 ~ —p^ ~ —V(<p). 
In this environment we may, to first approximation, ignore the last three terms in the right- 
hand side of f)7.5.ip . which then linearizes to 



A a = -6HA a , (7.5.6) 

where H = a/ a ~ constant and a cx exp(Ht). This result suggests an exponential decay for 
all three types of density inhomogeneities irrespective of the magnetic presence and in line 
with the cosmic no-hair theorems |161| - [T64"] . 

We note that here we have only considered the dynamical effect of the exponential expansion 
on the evolution of perturbations in the presence of the 5-field. Changes in the electrical 
properties of the cosmic medium, which can occur during the de Sitter regime, have not been 
accounted for. Also, typical linear studies, in magnetic-free universes, consider the second- 
time derivative of the density gradients instead of the first. This means that one decouples 
the expansion gradients from (17.5. ip before setting p + p ~ 0. One could do the same in the 
magnetic case as well. However, the calculation is rather involved and goes beyond the scope 
of this section. We expect that the inhomogeneities will still be found to decay exponentially, 
though probably at a slightly slower rate. 



8 Isocurvature perturbations 



One can define as isocurvature perturbations those occurring on hypersurfaces of uniform 
curvature, namely fluctuations which maintain D a 7l = at all times |147j . This should be 
distinguished from the definition typically found in the literature, where the term isocurvature 
means distortions in mult i- component systems with zero perturbation in the total energy- 
density initially (e.g. see |lb'5lll67j ). Here we will impose the T) a TZ = condition throughout 
the evolution of the perturbed mode, noting that one could use the equations given in section 
§ [7] to study fluctuations with an isocurvature initial condition by setting A + (A)/2)£> = 
(see §03). 



8.1 Consistency condition 



Isocurvature perturbations also require zero vorticity to guarantee the integrability of the 3-D 
hypersurfaces orthogonal to u a . In the magnetic presence, linear vortices are switched off by 
imposing the (self-consistent) constraint B b D b cvj:\B a = (see Eq. (16.2.41) ). The zero- curvature 
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condition for isocurvature magnetic perturbations is obtained through the orthogonally pro- 
jected gradient of (I6.2.9p . To be precise, using definitions (I5.7.ip . (15.7.2a ) and linearising 
around our FRW background we arrive at 

aD a K = 2pA a + clp(l + w)B a - AHZ a . (8.1.1) 

When D a 7l = 0, the right-hand side of the above ensures that linear expansion gradients 
and those in the fluid and the magnetic energy densities are connected by a simple algebraic 
relation. The projected comoving divergence of the latter translates into the following linear 
constraint between the associated scalar variables 



2HZ = p 



[8.1.2) 



In the absence of the magnetic field, this condition is automatically satisfied for pressure-free 
dust , but holds on large scales only when the matter has non-zero pressure [147] . If a magnetic 
field is present the vanishing of the matter pressure is not enough to guarantee that T) a TZ = 
at all times, because of the Lorentz-force contribution to the 4-acceleration. In that case the 
consistency condition for linear isocurvature perturbations is satisfied on large scales only. 
Indeed, on using the linearized propagation equation of the spatial Ricci scalar, we arrive at 

(Pan)' = -3 ( 1 + t cl) HDJZ + -4^- (l - ~ cl) D 2 A a + ^ D 2 B a , (8.1.3) 



9 V " a(l + w)\ 3 V " a 

since the vorticity has already been switched off. This result shows that the linear isocurvature 
condition is self-maintained only asymptotically (i.e. at infinity) where the Laplacians of the 
right-hand side vanish. Nevertheless, following |147| - [14"9"] . we will assume that on sufficiently 
long wavelengths the source terms in the right-hand side of (18.1.31) are negligible. Note the 
last term in the right-hand side of the above, which demonstrates why zero fluid pressure does 
not automatically guarantee the consistency of D a n — in the presence of the 5-field. 

Using the isocurvature condition (18.1.21) . we can eliminate Z from Eq. (17.2.11) . which reduces 
to 

A = -1 (l - w + - c>) HA + \ c 2 a (l - w*)HB - cl(l + w)HK . (8.1.4) 

The above, together with (17.2.31) and (17.2.4)) . describes the linear evolution of isocurvature 
scalar perturbations on a weakly magnetic flat FRW background filled with a single highly 
conducting perfect fluid. 



8.2 Evolution in the radiation era 



Since isocurvature perturbations apply to super-horizon scales only, all gradients in the fluid 
and the magnetic pressure have been switched off. This means that the only support against the 
gravitational pull of the matter comes from the expansion of the universe. When relativistic 
matter dominates the energy density of the latter, the system (I8.1.4p . (17.2.3)) and (17.2.41) 
reduces to 
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A= ~l{ 1 + l p ) r ' A ~ \ p r%K + 1 pt ' 1B ' (8 - 2 - 1} 

t = -~ /3r x /c + i (1 - ~ /?) r x A + 5 /^s , (8.2.2) 

B = A. (8.2.3) 

Taking into account that /3 = constant <C 1 and keeping up to /3-order terms, the solution for 
magnetised isocurvature perturbations reads 

A = 1 /3(A - So) + Ci^ + Cat-t 1 -^/ 2 , (8.2.4) 

where A , B are the perturbed matter and magnetic energy densities initially and Ci j2 are 
constants (see also [101J). When A ^ B , the isocurvature density contrast approaches a 
non-zero value which depends on the initial conditions. 



8.3 Evolution in the dust era 



After matter-radiation equality, the system (18.1.41) . (17.2.31) and (17.2. 4p becomes 



A = -r 1 a - 1 /^/c + 1 ^r 1 ^ 

9 3 



1a. 

3 



(8.3.1) 
(8.3.2) 
(8.3.3) 



with (3 = Po(to/t) 2 ^ 3 throughout the dust era. Because of this time- variation in the effective 
sound speed, the above system has no simple analytical solution. Nevertheless, taking the time 
derivative of (18.3. If) . using (18.3.21) . (18.3.31) and keeping up to /?o-order terms we arrive at 



A = — 



t~ 2 A. 



5.3.4) 



At late times, we may ignore the weak and decaying /3o(to/t) 2 ^ 3 term in the right-hand side of 
the above, which then admits the power-law solution [101J 



A = &r 2/3 + c 2 r x . 



(8.3.5) 



Comparing results (18.2.41) and (18.3.51) to their magnetic-free counterparts, shows that the 
presence of the -B-fieldd has only added new decaying isocurvature modes. Recall that in 
single-fluid almost-FRW cosmologies, isocurvature modes evolve as A oc t~ 1 ^ 2 during the 
radiation era and like A oc t^ 1 after equality [147J. 
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9 Entropy perturbations 



9.1 Entropy perturbations in multi- component systems 

Entropy perturbations can and will generically arise in all multi-component systems. The mag- 
netized, highly conductive barotropic medium may be treated as a mixture of two comoving 
fluids, of which one (the matter component) is perfect and the other (the I?-field) is imperfect 
(see Eq. (15.3.11) ). Then, it would be interesting to see whether some of the magnetic effects 
can be interpreted as coming from an effective entropy perturbation caused by differences in 
the dynamical behaviour of the two fluids. 

Assuming that £® is the intrinsic entropy perturbation associated with the i-th component 
of a multi-fluid system, the effective entropy perturbation of the total fluid is given by the 
expression [125] - |127] 

£ a = -Ep^P + ^tE h®hW ( c f ) - <*>>) . (9.1.1) 

Here, h® = p® + p( l > represents the inertial mass of the species, while p = Ya P an d 
h = J2i ^ are respectively the isotropic pressure and the inertial mass of the mixture. Also, 
the orthogonally projected vector 

c{ij) _ a (i) _ Pf!_ a (j) (q i o) 

with A™ = (a/p( l >)D a pW and = —S^ l \ describes the effective entropy perturbation 
triggered by the different dynamical behaviour of the individual components. Therefore, in 
a multi-fluid system, the total entropy fluctuation has one part coming from the intrinsic 
entropy perturbation of each species and one coming from differences in their dynamics. In 
this sense, a state of effective overall adiabaticity corresponds to £ a = 0. Note that, according 
to (19.1.11) . there is no entropy perturbation coming from differences in the dynamics of the 
member species when these share the same effective sound speed. 

9.2 Magnetically-induced entropy perturbations 

Treating the highly conducting magnetized medium as a system of two comoving fluids, we 
may only consider entropy perturbations due to the different dynamical evolution of the two 
components. These are determined by 

c(i2) = 1 A --B (9 2 1) 

given that p« = p, p^ = B 2 /2, /i« = p + p, = 2B 2 /3, c 2 ^ = c 2 s , c 2 ^ = 1/3 (since 
p B = pB /3 = £2/6), A« = A Q , =B a ,p = p+B 2 /Q and h = p+p+2B 2 /3. However, linear 
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perturbations in the magnetic energy density evolve in step with their matter counterparts 
(see Eq. (16.3.191) ). This guarantees that 



S£ 2) = \c a , (9.2.2) 

where C a = B® — 4A°/3(l + w) is a constant that depends on the initial conditions. Substituted 
into definition (19.1.11) . the above gives 

3p(l + W )(3cg-l)i? 2 
a ~ (6p + B 2 ) [3p(l + w) + 2B 2 } a ' 1 } 

This vanishes during the radiation era - a result guaranteed by the radiation-like behaviour 
of the .B-field (i.e. by the fact that p B = Pb/3 - see also [166] ). In the dust epoch, on the 
other hand, we find that £ a ~ C a = constant (since B 2 <C p always). These mean that, 
within the two-fluid description, part of the magnetic effect on the linear evolution of matter 
inhomogeneities can be interpreted as an effective entropy perturbation, due to the different 
dynamics of the two perturbed fluids, but only after equality. Then, the value of the effective 
entropy fluctuation depends on the initial relation between the magnetic and the fluid energy- 
density perturbations. In particular, when B a = 4A a /3 initially the dynamically induced 
entropy perturbation vanishes and we have effective adiabaticity. 



10 Vector perturbations 



A general inhomogeneous perturbation is characterized by its rotational and deformation 
properties, in addition to the amount of matter aggregation (or dilution). Following [lUlj . we 
will consider here the effect of magnetic fields on the linear evolution of vortex-like distortions 
in the density distribution of the cosmic medium. 



10.1 Basic variables 



Rotational instabilities in the fluid density are described by means of the antisymmetric, 
comoving, orthogonally projected tensor W a b = aD[&A a ] (see decomposition (16.3. 7|) in § !6.3.3p . 
The antisymmetry of the latter means that we can use it to define the vector 

W a = l -e abc W bc = -^curlA a , (10.1.1) 

with D a W a = to linear order. Recalling that A a = (a/p)D a p and using the general relativistic 
commutation law D[ a D&]</> = —4>ui a b between the projected gradients of scalars (see Eq. ( 1A.1I) 
in §E3J, we find that 

W a = -3a 2 (l + w)Huj a . (10.1.2) 
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This result reflects a fundamental property of vortex-like distortions in general relativity, 
namely that they are proportional to the amount of vorticity. Also, (I10.1.2P and the linearized 
version of constraint (12.3.61) guarantee that W a is a solenoidal vector to first order. Similarly, 
linear vortices in the volume expansion and the magnetic energy density are also directly 
related to the vorticity vector and therefore to those in the density of the matter. In particular, 

^curlZ Q = W a and %urlff a = - — 1 w , (10.1.3) 

2 (1 + w)H 2 3(1 + w) 

to first approximation. Finally, combining the antisymmetric part of the linearized Gauss- 
Codacci equation (see (15.6.31) in § 15. 6p with Eq. (110. 1.2j) . we obtain 

n a = - 1 w a , (10.1.4) 

3a 2 (l + w) 

where lZ a = e a b c 7l bc /2 is the vector component of the local spatial curvature. 



10.2 Linear equations 



Magnetic fields are unique in the sense that they are the only large-scale vector source that has 
ever been observed in the universe. Given that, the study of vortex-like, vector perturbations 
in the presence of the 5-field acquires particular significance. In our analysis, effects directly 
related to the vector nature of magnetic fields propagate through the tension part of the 
Lorentz force (see decomposition (15.3.41) ). Thus, as far as the scalar/density perturbations 
are concerned, these effects are encoded in the magneto-curvature terms seen in Eqs. (I7.2.ip . 
(J7X2D. 

Spatial inhomogeneities in the density distribution of a magnetized and highly conductive 
barotropic perfect fluid propagate according to the linear expression (16.3.111) in § 16.3.41 When 
combined with (16.3.141) . the latter takes the form 



A a = 3w (l - H A HA a - (1 + w)Z a + ^ cl(l + w)HB a 

-^(l- 2 -ci)B b D b B a , (10.2.1) 

where we have decomposed the magnetic Lorentz force into its isotropic-pressure and tension 
parts (see (15.3.41) in § 15. 3j) and kept up to c^-order terms. Taking the comoving, orthogonally 
projected gradient of the above and isolating its antisymmetric part we arrive at 



2p 



3 Ca 



[l-3wY-12w \HW a 



B b D b cm\ B a 



(10.2.2) 
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to linear order. Note that in deriving this expression we have also used the auxiliary relations 
fU0.1.3p . fll0.1.4p . the commutation laws between the spatial gradients of scalars and orthog- 
onally projected vectors and the zero-order versions of Eqs. (I5.5.ip . (I5.6.4p . The above shows 
that the -B-field acts as a source of density vortices when -B fe D&curl-B a ^ 0. Not surprisingly, 
the same requirement can also lead to kinematic vorticity (see Eq. (16.2.41) in § I6.2p . In fact, 
recalling that w = c 2 = constant, one can recover (16.2.41) from (110.2.21) by means of relation 
(I10.1.2p . An additional magnetic effect is the precession of W a , since its no longer parallel to 
W a (see also [TUT] ). 

For comparison reasons it will help to consider first the solution of Eq. (110. 2. 2p in the absence 
of the .B-field. Assuming conventional matter and dropping the magnetic terms from the right- 
hand side of the above, it is straightforward to show that non-magnetized fluid vortices decay 
always and on all scales. In particular, W a oc t~ 1 / 2 during the radiation era and ~W a oc t^ 1 
after equality. In what follows we will illustrate how the magnetic field changes this picture 
by focusing on the dust era, referring the reader to ["101 J for further discussion and details. 



10.3 Evolution in the dust era 

Expression (110.2.21) shows that the presence of the -B-field will generally trigger vortices in the 
density distribution of the magnetized matter. Nevertheless, the magnetic effect on preexisting 
rotational distortions is not clear yet. To quantify the role of the field we need to go one step 
further and obtain a decoupled equation for the evolution of W a . After equality, the Alfven 
speed decays with time according to c 2 = (3 oc t~ 2 ^ 3 . This means that, very soon, the weak 
Alfven terms in the right-hand side of (110.2.21) become completely negligible. At this limit, 
the time derivative of the latter leads to [101] 

W a = -AHW a - l - pW a + ± c 2 J) 2 W a . (10.3.1) 

Thus, during the dust era, magnetized vortices propagate like Alfven waves with signal speed 
Va, = c a /v3. Using standard vector harmonics, we may decompose the solenoidal W a as 
>V a = J2kW{k)Q { a k \ where D a W (fc) = = Q« = D a g« and D 2 Q( fc ) = -(k/a) 2 Q^. Then, 
recalling that w = = c 2 , a oc t 2//3 , H = 2/3t and p = 4/3t 2 , the k-th. Fourier mode propagates 
according to 

yy« = _^-iyv« - - 
3 3 

with Afc = a/k representing the wavelength of the perturbation and A a = ^/]3Xh- This repre- 
sents a characteristic scale of the magnetized environment, which may be termed the 'Alfven 
horizon' [101] . The above admits the solution 

W (fe) = dt ai + C 2 t a2 , (10.3.3) 



1 + 
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t -2 W (fc) 



(10.3.2) 
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where Ci,2 are constants and 
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(10.3.4) 



On scales far beyond the Alfven horizon, X a /Xk *C 1 and the perturbed mode exhibits a 
power-law decay with W oc t~ 2 / 3 , a rate considerably slower than the one associated with dust 
dominated, magnetic-free cosmologies. Recall that W oc t" 1 in the absence of the 5-field. Well 
inside A a , on the other hand, perturbations oscillate like Alfven waves with a time-decreasing 
amplitude. In particular, for A a /Afc ^> 1, solution (110.3.31) . (I10.3.4p implies 



W(fc) oc cos 



(±) In* 



-5/6 



(10.3.5) 



Thus, the overall, the effect of the I?-field on a given vortex mode is to decrease its standard 
depletion rate. 



Not surprisingly, a similar magnetic effect is also observed on vorticity proper. Using the 
linear relation ( 110.1. 2ft . we find that after equality large-scale rotational distortions decay as 
u a oc t _1 as opposed to uj a oc t~ 4//3 . On these grounds, magnetized cosmologies would contain 
more residual vortices than their magnetic-free counterparts. 



11 Gravitational waves 



Gravitational waves are propagating fluctuations in the geometry of the spacetime fabric, 
usually described as weak perturbations of the background metric. Alternatively, one can 
describe gravity waves covariantly by means of the electric and magnetic components of the 
Weyl tensor |168j . which encodes the long-range gravitational field, namely the one determined 
by the presence of matter "elsewhere" in the space (see § 14.21) . 



11.1 Isolating magnetised tensor modes 



Covariantly, gravitational waves are described by the transverse degrees of freedom in the 
electric (E ab ) and magnetic (H ab ) parts of the conformal curvature tensor. The transversal- 
ity is necessary to ensure that the pure tensor modes of the locally free gravitational field 
have been isolated. The same condition is also imposed on the shear and any other orthogo- 
nally projected, traceless, second-rank tensor that might be present. Thus, when studying the 
propagation of gravitational radiation in perturbed FRW models with perfect fluid matter we 
require that (see § 12.31 and § 14.21) 

D b E ab = = D b H ab = D b a ab , (11.1.1) 
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to linear order and at all times (e.g. see [169j ). In magnetic-free universes, this is achieved by 
switching the vorticity off and by setting D a p = = D a p = D a B (for a barotropic medium it 
suffices to ensure that D a p = = D a 0). These constraints, which are self-consistent (i.e. pre- 
served in time) at the linear perturbative level, guarantee that the 4-acceleration also vanishes 
to first approximation. When a magnetic field is included there is an additional constraint on 
the anisotropic pressure of the field. In particular, the system (jll.l.ip reads [14211170] 



D b E ab = \ D aP + \ D a B 2 - I D b U ab = , (11.1.2) 
6 b 2 

D b H ab = p(l + w)u a = 0, (11.1.3) 

D 6 a ab = ^D a e + curl^ a = 0, (11.1.4) 

D b U ab = e abc B b cm\B c - \ D a B 2 = . (11.1.5) 

o 



When studying the magnetic effects on cosmological gravitational waves, one needs to take 
into account the anisotropic nature of the 5-field. It is therefore not appropriate to assume 
a fully random 5-field, which simply acts as an additional source of energy density. Thus, 
following [68] and also [14211170] , we will treat both the energy density and the anisotropic 



pressure of the field as first order perturbations|__J This approach differs from the perturbative 



scheme adopted for the scalar and the vector modes, where there was a weak 5-field in 
the background. Nevertheless, adopting the scheme of [SB] simplifies the mathematics of the 
gravito-magnetic interaction without compromising the physics. 

We assume a spatially flat FRW background. Then one can isolate the pure-tensor pertur- 
bations by adding to the previously mentioned perfect-fluid constraints the following linear 
conditions [142lf l70j 

B a B 2 = = e abc B b cm\ B c . (11.1.6) 

In other words, both magnetic energy-density gradients and the Lorentz force vanish at all 
times. These guarantee that there is no magnetic contribution to the 4-acceleration and sub- 
sequently that the field does not trigger any vorticity, density or expansion perturbations. The 
consistency of constraints (111.1. 6p is straightforward to show. In particular, we find that 

(D a B 2 y = -5H(D a B 2 ) and (e abc B b cm\ B c \ = -5He abc B b cm\ B c , (11.1.7) 

to lowest order. This means that the constraints (111.1.61) continue to hold after they are 
initially imposed. Then, the only remaining nontrivial linear constraints are 

H ab = curl a ab and TZ {ab) = -Ha ab + E ab + - U ab , (11.1.8) 



13 F orm ally, treating B 2 = B a B a as a first-order variable, means that B a is of order 1/2 in perturba- 
tive terms. Then, because T) a B 2 is of order 1 and D a l? 2 = 2B b D a B b , consistency demands that the 
projected magnetic gradients are also half-order perturbations. 
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where curl <?ab = e cd{a D c a b) d (see Eqs. (12XTD and (14X61) in §Pand § Irrespectively). Note 
that, according to (111. 1.8b ). the linear conditions D b E ab = = D b a ab = D b H a b guarantee that 
D b TZ {ab) = as well. 



11.2 Covariant gravitational-wave energy density 



The energy density of gravitational radiation is determined by the pure tensor, namely the 

TT- 

x/3 . 



transverse traceless part (h^T) of the metric perturbation, according to (e.g. see |171j ) 



Pgw = 2^ , (11.2.1) 

where a prime indicates differentiation with respect to conformal time (recall that c = 1 = 87rG 
throughout this review). In a comoving frame, with u a = 5o a u°, we have a 2 = a ab o ah jl = 
cr Q/ 3<7 a/3 /2, with the transverse part of the shear components given by 



cr. 



ap = a{ti%)' and = a" 3 (/^)' . (11.2.2) 

Solving these relations for the metric perturbations and then substituting the results into 
Eq. flll.a.ip . we arrive at [H2] 

Pgw = o- 2 } (11.2.3) 

which provides a simple covariant expression for the energy density of gravitational-wave 
distortions. 



11.3 Evolution equations 



In an FRW spacetime the Weyl tensor vanishes identically, which means that E ab and H ab 
provide a covariant and gauge-invariant description of perturbations in the free-gravitational 
field. Once the pure tensor modes have been isolated, we can proceed to linearize the prop- 
agation equations (I4.2.5p . (I4.2.6P of § 14.21 On a Friedmann background with a single perfect 
fluid, the latter reduce to 



E ab = -3HE ab - - p(l + w)a ab + cm\H ab - ift^ - 1 HU ab , (11.3.1) 
H ab = -3HH ab - cm\E ab + l - curlIT afe . (11.3.2) 

Since the magnetic part of the Weyl tensor satisfies the constraint (111.1.8a ). the linear evolu- 
tion of H ab is determined by the shear evolution (12.3.31) 



& ab = -2Ha ab -E ab + ^U ab . (11.3.3) 
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Furthermore, on using the commutation law between the orthogonally projected gradients of 
spacelike tensors, constraint (111. 1.8b ) leads to the auxiliary relation curlif a f, = — D 2 a ab , and 
Eq. (111.3. ip becomes 

E ab = -3HE ab - - p(l + w)a ab - ^tl ab - ^ HU ab - D 2 a ab . (11.3.4) 

This, together with Eq. (111.3.31) and the linearized part of (15.2.21) . namely 

tl ab = -4ittl a6 , (11.3.5) 

describe the linear evolution of gravitational waves in almost-FRW universes containing a 
single highly conductive barotropic fluid and a magnetic field. 

We can describe this evolution in terms of the shear as follows. Taking the time derivative 
of (111.3.31) . using Eqs. (jll.3.4p . (I11.3.5p . the background Raychaudhuri and Friedmann equa- 
tions, and keeping only linear order terms, we arrive at the following wave equation for the 
gravitationally induced shear |170] 

a ab - D 2 cr a6 = -5H& ab - i p(l - 3w)a ab - 2HIl ab . (11.3.6) 

This is no longer coupled to the electric Weyl tensor, and together with Eq. ( 111. 3. 5ft . describes 
the propagation of gravitational waves. The magnetic anisotropic stress can act as a source for 
gravitational waves, and it is possible to place strong constraints on magnetic fields via this 
generation of gravitational waves |171lll73j . Also, the anisotropic stress of the B-field slows 
the decay of shear |174] . 

Following (lll.3.6p . the solution of Eqs. (111.3.51) and (I11.3.6P depends critically on the equation 
of state of the matter component, and in particular on whether w = 1/3 or not. In the matter 
dominated era and on superhorizon scales, we find 

a = p 1 r 1/3 + p 2 r 2 + Qr 5/3 , (11.3.7) 

where Pi = = P2 = Q. The P-modes are the usual Bianchi / type solutions, while the 
Q-mode carries the effect of the I?-field. The latter mode means that, although the magnetic 
presence does not alter the standard evolution rate of the shear, it can affect its magnitude 
(in a rather involved way - see [142J). Note that in deriving the above we have also introduced 
the familiar tensor harmonics Q^ b \ with = QS)> Q^ab = = D fe Q^ and D 2 Q^ = 
-{k/a) 2 Q { ^. These were used to decompose a ab and Ii ab according to a ab = J2k a {k) Q^ab an d 
n afe = E fc n (fc) Qi? respectively (with D a( x« = = D a Il( fc )). 

11.4 The zero- eigenvalue issue 

The evolution equation of trace-free anisotropic stresses during the radiation era presents a 
particular mathematical problem, when linearised around an FRW background, because of 
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the zero eigenvalue issue (see [65] and references there in). In particular, following Eq. (I6.1.6p . 
it appears that the Alfven speed and therefore the ratio B 2 / p remain constant as long as 
w = 1/3. This is a critical case, however, and a closer analysis leads to B 2 /p oc (logt) -1 
because of second-order pressure effects. When these are taken into account, we find that |174] 

(Tab = Pabt- 3/2 + Qabitlogt)- 1 , (11.4.1) 

on superhorizon scales. The magnetized mode falls off slowly due to the anisotropic pressures, 
with cr a 5 oc (tlogt) -1 , and there is an attractor with a/H — > B 2 / p and IT a b — > 2E a i, oc 
(t 2 logi) -1 when the dynamics of the radiation era are close to isotropy [65.174J. In this 
regime the leading-order diagonal scale factors of the spacetime metric during the radiation 
epoch take the form 

di(t) oct 1/2 (logt) 9 ', (11.4.2) 

with the {qi, i = 1, 2, 3} constants satisfying the constraint Z^ =1 q% = 0. In the dust-dominated 
era the anisotropies fall as power-law in time and the orthogonal scale factors evolve as 

ai (t) oct 2/3 (l + Vit- s >), (11.4.3) 

where the {Vi, i = 1,2,3} and the {s^, i = 1,2,3} are sets of constants with Y^=i s % — 0. 
Since the observed microwave background temperature anisotropy is determined by the value 
of a/H at recombination, and a/H oc B 2 / p m oc t~ 2 / 3 after the epoch of matter-radiation 
equality, we obtain a strong observational bound on the present-day magnetic field strength 
from the large-angular scale anisotropy of the microwave background at last scattering, when 
z ~ 1100, [66]: 

\B\ < 4.0 x 10 _9 fio /2 G. (11.4.4) 

An investigation of the behaviour of Yang-Mills fields in the early universe reveals that chaotic 
behaviour is possible in the evolution of the scale factors with time, no matter how weak the 
Yang-Mills field, and how close the expansion is to isotropy [175|I176] , if there is no perfect fluid 
present in addition to the magnetic field. However, despite the close relationship between the 
form of Yang-Mills fields and magnetic fields, the slow logarithmic decay of a/H oc (logt) -1 
of the shear distortion with time characteristic of magnetic fields and radiation does not occur 
if an isotropic blackbody radiation fluid is added to the Yang-Mills fields |177j . because of the 
different time-dependence of their pressures (Jl a b/PYM is no longer constant), which hastens 
their decay, and the shear falls off in Bianchi / universes at the same rate as it does when 
anisotropic stresses are absent, with a/H oc t~ l l 2 . This means that the microwave background 
anisotropy provides a rather weak bound on the energy- density in Yang-Mills fields today of 
Hyjv/0 < O.llllyo, where f2 7 o is the present density of blackbody fluids in units of the critical 
density. By contrast, the slow decay of the shear in the presence of magnetic stresses led to a 
far stronger bound on the magnetic field density today, with Q mag o < 10~ 5 fi 70 . 
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12 Summary 



Magnetic fields are detected everywhere when the appropriate observations are possible. The 
origin of cosmic magnetism is still an open issue but the presence of large-scale .B-fields with 
similar /zG-order strengths in galaxies, galaxy clusters, and high-redshift protogalactic struc- 
tures may suggest a common (primordial) origin for them. Large-scale, micro-Gauss fields can 
significantly affect the evolution of the universe and are extremely important for astrophysics 
(see [T2] for a cross-disciplinary review of the subject). The idea of primordial magnetism is 
appealing because it can potentially explain all the large-scale fields seen in the universe to- 
day, especially those found in remote protogalaxies. As a result, the literature contains many 
studies examining the role and the implications of magnetic fields for cosmology. 

In this work we have employed covariant techniques to analyse the effects of magnetic fields on 
the evolution of inhomogeneous relativistic cosmologies. After a brief introduction to nonlinear 
cosmological electrodynamics and magnetohydrodynamics, we investigated weakly magnetized 
almost-FRW universes. Within the standard ideal MHD limits, we found that linear inhomo- 
geneities in the magnetic energy density grow in step with those in the matter throughout the 
hot big-bang evolution of the universe. This close relation between the magnetic and the fluid 
density gradients is a key property of weakly-magnetized FRW cosmologies. 

We found that during the radiation era and on large scales the magnetic pressure inhibits the 
growth of density perturbations by an amount proportional to the field's relative strength. On 
sub-horizon lengths the extra magnetic pressure was shown to trigger magnetosonic waves with 
a reduced amplitude and an increased frequency, relative to the magnetic-free case. Moreover, 
the magnetic presence means that small-scale density perturbations no longer oscillate around 
a zero average value but have a finite residual average that depends on the initial conditions. 
After equality, the B- field is effectively the only source of pressure support. This reduces 
the growth rate of matter condensations in a way exactly analogous to that observed on 
super-Hubble scales during the radiation era. The same magnetic pressure can also stabilize 
against gravitational collapse by bringing the Jeans length up to the size of a galaxy cluster. 
Looking at models with a vacuum equation of state for the matter, we found that the outcome 
depends on the effective equation of state of the medium. In particular, the comoving density 
gradients remain constant as long as the barotropic p = —p condition holds. When dealing 
with non-barotropic minimally coupled scalar fields, on the other hand, all three types of 
density inhomogeneities suffer exponential decay. 

Defining as isocurvature large-scale perturbations that maintain the zero 3-curvature require- 
ment at all times, we found that they contain decaying modes only. Given that, typically, the 
term 'isocurvature' refers to fluctuations with zero initial total energy- density perturbation, 
we have explained how our equations can be used to study this type of perturbation mode. By 
treating the highly conductive medium as a two-component system we found that the effective 
entropy perturbation caused by the different dynamical evolution of the two species is either 
constant or zero. The evolution of vector perturbations in the presence of magnetic fields, 
showed that the latter increase the amount of residual vorticity. Thus, universes containing 
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large-scale inhomogeneous 5-fields would rotate more than their magnetic-free counterparts. 
Finally, we investigated the implications of the field for the evolution of gravitational-wave 
perturbations. 
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A Commutation laws 

According to definition (12.1.2b ,) . the orthogonally projected covariant derivative operator sat- 
isfies the condition D a hb c = 0. This means that we can use h a b to raise and lower indices in 
equations acted upon by the aforementioned operator. Following Frobenius' theorem, however, 
rotating spaces do not possess integrable 3-D submanifolds (e.g. see [14011141] ) . Therefore, the 
D a -operator cannot be used as a standard 3-D derivative in such spaces and it does not always 
satisfy the usual commutation laws (see below and also |149j ). 



A.l Orthogonally projected gradients 

When acting on a scalar quantity the orthogonally projected covariant derivative operators 
commute according to 

D [a D 6] ./= -Uabf. (A.l) 

The above is a purely relativistic result and underlines the different behaviour of rotating 
spacetimes within Einstein's theory. Similarly, the commutation law for the orthogonally pro- 
jected derivatives of spacelike vectors reads 

D [a D b] v c = -uj ab v {c) + - TldcbaV d ■ (A. 2) 

where v a u a = and TZ a bcd represents the Riemann tensor of the observer's local rest-space. 
Finally, when dealing with orthogonally projected tensors, we have 

1 

DfoDbJiS'cd = —^abhchd S e f + — (lZ e cbaS e d + T^edbaSc) , (A. 3) 

with S a bU a = = S a bU b . Note that in the absence of rotation, TZ a b c d is the Riemann tensor of the 
(integrable) 3-D hypersurfaces orthogonal to the w a -congruence For details on the definition, 
the symmetries and the key equations involving TZ a b c d, the reader is referred to § 14.31 We also 
note that the above equations are fully nonlinear and hold at all perturbative levels. 
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A. 2 Time derivatives 



In general relativity, time derivatives do not generally commute with their spacelike counter- 
parts. For scalars, in particular, we have 

DJ - h a \D b fy = -fA a + 1 0D a / + D b f [a\ + u\) , (A.l) 

at all perturbative levels. Assuming an FRW background, we find that the orthogonally pro- 
jected gradient and the time derivative of the first-order vector v a commute as 

aD a v b = (aD a v b )' , (A. 2) 

to linear order. Similarly, when dealing with first-order spacelike tensors, we have the following 
linear commutation law 

aD a S bc = (aD a S bc )' . (A. 3) 



B Notation 

• Spacetime Geometry 

Line element: ds 2 = g ab dx a dx b = — dr 2 , with c = 1. 
4- velocity: u a = drr a /dr, 3-D projection tensor: h ab = g ab + u a u b . 
4-D permutation tensor: r] abcd , 3-D permutation tensor: e abc = i] abcd u d . 
Covariant derivative: V b T a = dT a /dx b — T ab T c . 

Time derivative: T a = u b, V b T a , 3-D covariant derivative: D b T a = h b d h a c V dT c . 

Riemann tensor: R a bcd, Ricci tensor: R ab = R c ac b, Ricci scalar: R = R a a . 

3-Riemann tensor: 7Z abc d, 3-Ricci tensor: 1Z ab = 1Z c acb) 3-Ricci scalar: 1Z = lZ a a . 

3-curvature index: K = 0, ±1, with TZ = 6K/a 2 (in FRW models). 

Weyl Tensor: C abcd , electric Weyl: E ab = C acbd u c u d , magnetic Weyl: H ab = e a cd C cdbe u e /2. 

• Kinematics 

Expansion scalar: = V a- u a = D a u a , scale factor: a, with a/a = 0/3. 
Vorticity tensor: u ab = D^u a ], vorticity vector: u a = e abc u bc /2. 
Shear tensor: a ab = D^w^ = D( b u b ) — (D c u c )h ab /3, 4-acceleration: A a = u b V b u a . 
Hubble parameter: H = a/a (in FRW models). 

• Matter Fields 

Field equations: R ab — (R/2)g ab = T ab , with k = 8irG = 1. 

Matter energy-momentum tensor: T ab = pu a u b + ph ab + 2u( a q b ) + n ab . 

Matter density: p = T ab u a u b , isotropic pressure: p = T ab h ab /3. 

Barotropic index: w = p/p, adiabatic sound speed: c 2 = dp/ dp. 

Energy flux: q a = h a b T bc u c , anisotropic pressure: n ab = T( a6 > = T {ab ) - (T/3)h ab . 

• Electromagnetism 

Electromagnetic tensor: F ab , magnetic field: B a = e abc F bc /2, electric field: E a = F ab u b . 
Magnetic energy density: B 2 /2, magnetic isotropic pressure: B 2 /Q. 
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Magnetic anisotropic pressure: U a b = —B( a Bq = (B 2 /3)h a b — B a Bb. 
Alfven speed: c 2 & = B 2 /(p + p + B 2 ). 

Electric 4-current: J a , electric 3-current: J a = J( a ) = h a b Jb- 
Charge density: p e = —J a u a , electrical conductivity: q. 
• Perturbations 

Matter density gradients: A a = (a/p)D a p, with A a b = aD b A a and A = A a a . 
Matter vortices: W ab = A [ab] , with W a = e abc W hc /2. 

Volume expansion gradients: Z a = aD a 0, with Z ab = aD b Z a and Z = Z a a . 
Magnetic density gradients: B a = (a/B 2 )D b B 2 and B a b = aD b B a , with B = B a a . 
Effective entropy perturbations: £ a , with 5™ = —S^ % \ 
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